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PREDICTION OF VORTEX SHEDDING FROM 
CIRCULAR AND NONCIRCULAR BODIES IN 
SUBSONIC FLOW 

i 

Michael R. Mendenhall 
Daniel J. Lesieutre 

i 

I 

Nielsen Engineering & Research, Inc, 

| SUMMARY 

| An engineering prediction method and associated computer 

code VTXCLD to predict nose vortex shedding from circular and 

1 

i noncircular bodies in subsonic flow at angles of attack and roll 

are presented. The ax i symmetric body is represented by point 
i sources and doublets, and noncircular cross sections are trans- 

formed to a circle by either analytical or numerical conformal 
transformations. The lee side vortex wake is modeled by discrete 
vortices in crossflow planes along the body; thus the three- 
dimensional steady flow problem is reduced to a two-dimensional, 
unsteady, separated flow problem for solution. Comparison of 
measured and predicted surface pressure distributions, flow field 
surveys, and aerodynamic characteristics are presented for bodies 
with circular and noncircular cross sectional shapes. 

INTRODUCTION 

Current missile applications reguiring high aerodynamic 
performance can involve noncircular body shapes in subsonic flow 
at high angles of attack and nonzero roll angles. The angle of 
attack range may be sufficiently high to cause formation of body 
separation vortices, and the body vortex shedding characteristics 
are directly influenced by the body cross-sectional shape and the 
flow conditions. It is desirable to model the body vortex wake 
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by means of a rational method capable of considering a variety of 
body shapes over a wide range of incidence angles and Mach num- 
bers up to the critical speed. It is important that the separa- 
tion vortex wake induced effects on the nonlinear aerodynamic 
characteristics of the body be handled properly with a method 
which correctly represents the physical characteristics of the 
flow field. 

The phenomena of interest are the sheets of vorticity formed 
on the lee side of the body at high angles of attack. The vorti- 
city is formed by boundary-layer fluid leaving the body surface 
from separation points on both sides of the body (Fig. 1). At 
moderate angles of incidence, the vorticity rolls up into a sym- 
metrical vortex pair, but at higher angles, the vorticity becomes 
asymmetric as shown in the sketch in Figure 1(b). A method to 
predict these flow phenomena in the vicinity of circular and 
noncircular bodies in subsonic flow is described in References 1 
and 2. The extension of the method to supersonic flow is de- 
scribed in Reference 3. 

The purpose of this report is to describe an engineering 
prediction method and associated computer code developed to cal- 
culate the nonlinear aerodynamic characteristics and flow fields 
of noncircular bodies at high angles of attack at speeds up to 
the critical speed. The objectives of the method are to use a 
three-dimensional, attached flow, potential method to represent 
the body and a two-dimensional, incompressible, separated flow 
model to calculate the lee side vortex shedding from the body 
alone at angle of attack and angle of roll. The predicted pres- 
sure distribution on the body under the influence of the free 
stream and the separation vortex wake is used to calculate the 
aerodynamic loads on the body. Conformal mapping technigues are 
used to transform noncircular cross sections to a circle for 
calculation purposes. 
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The following sections of this report include a discussion 
of the approach to the problem and a description of the analysis 
and flow models required to carry out the calculation. The pre- 
diction method is evaluated through comparison of measured and 
predicted results on a variety of body shapes, including circular 
and elliptical cross sections. A user's manual for the computer 
code is also included. The manual consists of descriptions of 
input and output and of sample cases. 
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LIST OF SYMBOLS 
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C 


Pi 
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N 
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semi-major or horizontal axis of elliptic cross 
section 

coefficients of conformal transformation, Eq. (8) 

semi-minor or vertical axis of elliptic cross sec- 
tion 

normal-force coefficient per unit length, Eq. (49) 
side-force coefficient per unit length, Eq. (53) 
axial-force coefficient, Eqs. (58) and (62) 
rolling-moment coefficient, Eq. (57) 
pitching-moment coefficient, Eq. (51) 
yawing-moment coefficient, Eq. (55) 
pressure coefficient, Eq. (39) 
incompressible pressure coefficient, Eq. (41) 
side-force coefficient, Eq. (55) 
normal-force coefficient, Eq. (50) 

diameter of circular body, or equivalent diameter 
of noncircular body 

complex velocity component. Eg. (28) 
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*ref 

L 

M x 

M y 

M z 

M 

oc 

p 

Poc 

Pec 

0 

r 

r' 


R e 

S 


reference length 

total number of Fourier coefficients used to de- 
scribe transformation, Eg. (8); also model length 

rolling moment about the x-axis 

pitching moment about the y-axis 

yawing moment about the z-axis 

free-stream Mach number 

local static pressure 

free-stream static pressure 

1 2 

free-stream dynamic pressure, ^ pV w 
source strength 

radial distance from a vortex to a field point 

radial distance to a point on a noncircular body. 
Fig. 4 

vortex core radius, Eq. (59) 
radius of circle 

Reynolds number based on boundary layer run length 

and minimum pressure conditions, U £/v 

m 

reference area 
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V R 
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V ,w 


V ( 
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x,y,z 


x 


m 


a 


a 


c 


3 


3 1 


Y 


surface velocity in crossflow plane 
axial perturbation velocity 
local velocity 

radial perturbation velocity 
vortex induced velocity 
velocity components in real plane 
free-stream velocity 
complex potential. Eg. (26) 

body coordinate system with origin at the nose: x 

positive aft along the model axis, y positive to 
starboard, and z positive up 

axial location of center of moments 

angle of attack 

angle between free-stream velocity vector and body 
axis 

angle of sideslip; also polar angle in a-plane. 

Fig. 2; also „2 
1— M 

CO 

local slope of body surface. Fig. 4 
ratio of specific heats 
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Ax 
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n 

0 

v 
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P 
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T, X 


♦ 
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exterior angles of body segment for numerical map- 
ping, Eg. (6) 

vortex strength 

axial length increment 

complex coordinate in an intermediate plane. 

Fig. 2(b) 

vertical coordinate in an intermediate plane. 

Fig. 2(b) 

polar angle in v- plane. Fig. 2(a) 

complex coordinate in circle plane. Fig. 2(a); also 
kinematic viscosity 

run length, Egs. (46) and (47); or lateral coordi- 
nate in an intermediate plane. Fig. 2(b) 

free-stream density 

complex coordinate in real plane. Fig. 2(a) 

lateral and vertical coordinates in circle plane. 
Fig. 2(a) 

roll angle 

velocity potential in real plane 
stream function in real plane 
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Subscripts and Superscripts 


( ) conjugate of complex quantity 

( ) m vortex m 

( ) cp center of pressure 

( )' incompressible quantity; or surface values in 

Fig . 4 
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GENERAL APPROACH 


Bodies at high angles of attack exhibit distributed vorti- 
city fields on their lee side due to boundary-layer fluid leaving 
the body surface at separation lines. One approach for modeling 
these distributed vorticity fields has involved the use of clouds 
of discrete potential vortices. Underlying the basic approach is 
the analogy between two-dimensional unsteady flow past a body and 
the steady three-dimensional flow past an inclined body. In 
fact, the three-dimensional steady flow problem is reduced to the 
two-dimensional unsteady separated flow problem for solution. 
Linear theory for the attached flow model and slender body theory 
to represent the interactions of the vortices are combined to 
produce a nonlinear prediction method. The details of the appli- 
cation of this approach to the prediction of subsonic flow about 
circular and noncircular bodies are presented in References 1 and 
2. Other investigators have used this approach to successfully 
model the subsonic flow phenomena in the vicinity of circular 
cross section bodies (Refs. 4 and 5). The purpose of this report 
is to document in detail the subsonic analysis of References 1 
and 2 and to describe the code VTXCLD assembled to accomplish the 
subsonic vortex shedding calculations. 

The calculation procedure is carried out in the following 
manner. The axisymmetric body is represented by discrete point 
sources and doublets, and the strength of the individual singu- 
larities is determined to satisfy a flow tangency condition on 
the body in a nonseparated uniform flow at angle of incidence and 
roll. Compressibility effects on the body are included by a 
Gothert transformation which keeps the cross section shape un- 
changed but stretches the axial body coordinate. Starting at a 
crossflow plane near the body nose, the pressure distribution on 
the body is computed using the full compressible Bernoulli egua- 
tion. The boundary layer in the crossflow plane is examined for 
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separation using modified versions of Stratford's separation 
criteria. The Stratford separation criteria are based on two- 
dimensional incompressible flow. At the predicted separation 
points, incompressible vortices with their strength determined by 
the vorticity transport in the boundary layer are shed into the 
flow field. The trajectories of these free vortices between this 
crossflow plane and the next plane downstream are calculated by 
integration of the equations of motion of each vortex, including 
the influence of the free stream, the body, and other vortices. 
The pressure and trajectory calculations are carried out by map- 
ping the noncircular cross section shape to a circle using either 
analytical or numerical conformal transformations. The vortex 
induced velocity contribution to the body tangency boundary con- 
dition includes image vortices in the circle plane. At the next 
downstream crossflow plane, new vortices are shed, adding to the 
vortex cloud representing the wake on the lee side of the body. 
This procedure is carried out in a stepwise fashion over the 
entire length of the body. The details of the individual methods 
combined into the prediction method are described in the fol- 
lowing section. 


METHODS OF ANALYSIS 

The development of an engineering method to predict the 
pressure distributions on arbitrary missile bodies in subsonic 
flow at high angles of incidence requires the joining of several 
individual prediction methods. In this section, the individual 
methods are described briefly, and the section concludes with a 
description of the complete calculation procedure. 

Conformal Mapping 

The crossflow plane approach applied to arbitrary missile 
bodies results in a noncircular cross section shape in the pres- 
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ence of a uniform crossflow velocity and free vortices in each 
plane normal to the body axis. The procedure used to handle the 
noncircular shapes is to determine a conformal transformation for 
mapping every point on or outside the arbitrary body to a corre- 
sponding point on or outside a circular body. The two- 
dimensional potential flow solution around a circular shape in 
the presence of a uniform flow and external vortices is well 
known and has been documented numerous places in the literature 
(Refs. 6 and 7). Thus, the procedure is to obtain the potential 
solution for the circular body and transform it to the noncir- 
cular body. Conformal transformations used are of two distinct 
types, analytical and numerical. 

Analytic transformation .- For very simple shapes like an 
ellipse [Fig. 2(a)], the transformation to the circle can be 
carried out analytically as described in Reference 7. For 
example. 


where 


o 


v 



(1) 


= y + iz 


in the real plane, and 


v = x + i X 


(2) 


(3) 


in the circle plane. 


The derivative of the transformation is 



( 4 ) 
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which is required for the velocity transformation discussed 
later. 


Numerical transformation .- For complex noncircular shapes, 
the transformation cannot be carried out analytically and a nu- 
merical transformation is required. The numerical transformation 
chosen is described in detail in Reference 8, and a brief summary 
of the conformal mapping procedure follows. 

The sequence of events in the numerical mapping is shown in 
Figure 2(b). The arbitrary cross section shape of the body in 
the a-plane is required to have a vertical plane of symmetry. 
The transformation of interest will map the region on and outside 
the body in the a-plane to the region on and outside a circle in 
the v-plane. The first step is a rotation to the s-plane so 
that the cross section is symmetric about the real axis. Thus, 

C = i o 

= K + in (5) 

A mapping that transforms the outside of the body in 
the c-plane to the outside of the unit circle is 

. m 

' - I —2 ", (v 

v r=l 

where Y r are the exterior angles of the m segments of the body 
cross section. In Equation (6), n denotes a product series. For 
a closed body 

m 

l = 2tt (7) 

r=l r 


v ) 
r 


d v 


( 6 ) 
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The final transformation has the form 


a 


-i 




( 8 ) 


where the coefficients are obtained through an iterative 
scheme described in Reference 8 and r Q is the radius of the 
equivalent circle in the v-plane. The derivative of the trans- 
formation is 


do 
d v 


-l 


1 


L 

l 

1 = 1 


.*+1 


IA, 


£+1 


(9) 


which is required for the velocity calculations described in a 
later section. 

The above numerical mapping procedure has been applied to a 
wide range of general cross section shapes with good success. 

Body Model 

A three-dimensional representation of the missile volume is 
needed for purposes of predicting the absolute pressure coeffi- 
cient on the surface (Ref. 2). Since the model must be a compu- 
tationally efficient means of representing both circular and 
noncircular bodies in compressible flow up to the critical speed, 
a method using discrete singularities on the body axis is de- 
scribed in this section. As described in Reference 3, a panel 
method can be used to represent the missile surface; however, 
increased computational requirements make the present approach 
more desirable for an engineering prediction method. 

Circular bodies .- The volume of an axisymmetric body is well 
represented by a series of point sources and sinks distributed on 
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the axis. A number of models in varying degrees of complexity 
are available for this task; for example, see References 9-12. 
For use in VTXCLD, a discrete source model from the latter refer- 
ence was selected for its accuracy, economy, and reliabilty. 


Given a series of K point sources (sinks) distributed on the 
missile axis, the induced axial and radial velocities at a point 
( x , r ) are 



( 10 ) 


(ID 


where 


4irL V coso 
« c 

is the dimensionless source strength. Q£ and x^ are the source 
strength and axial locations, respectively, of the k-th source. 
The body surface slope at the j-th body point obtained from Egua- 
tions (10) and (11) is 


dr 

dx 


R . 


V cos a 


u 


1 + 


r . 


V COS a 


for j=l 


(N-3 ) 


(13) 
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Equation (13) produces a set of K-3 linear equations in the un- 
known source strengths. For a closed body, imposing the condi- 
tion of the sum of the source strengths to be zero gives the 
relation 

K 

I O' = 0 (14) 

k = l K 

The remaining two conditions are the imposition of stagnation 
points at the nose and tail of the body. These conditions from 
Equation (10) are 


K 

l 


k=l 



1 


(15) 



(16) 


Given the source positions on the missile axis. Equations (13) 
through (16) comprise the total set of eguations to solve for the 
source strengths. The predicted hull shape from the stream func- 
tion is 


* [f - 1 1 1 

r 

1 + ■ 

1 

X 1 " 
x k 

L “ L 


'[! - * 

[e] 2 " 

1/2 







which must be solved by iteration. 

The above method has proved successful in modeling a variety 
of axisymmetric missile bodies; however, some care is required 
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because of the ill-conditioned matrix. Best results are obtained 
if the discrete sources are spaced at intervals of 60% of the 
local radius and if the surface slope description of the body is 
smooth without discontinuities. Boundary conditions are satis- 
fied at points midway between the source locations. The predic- 
tion of an approriate source distribution is an automated part of 
VTXCLD, and the user must only specify missile geometry. 

Noncircular bodies .- An appropriate body model for missiles 
with noncircular cross sections is a surface panel method similar 
to that described in Reference 3; however, the use of a panel 
model adds significantly to the cost of each computation. For 
this reason an alternate approach was selected for use with non- 
circular bodies. 

The noncircular body is replaced with an eguivalent axisym- 
metric body having the same cross sectional area distribution as 
the actual body. There are approximations involved with this 
model since the induced u-velocity due to noncircular effects are 
neglected; however, for bodies that do not deviate greatly from 
circular, the approximations are acceptable. Some guidelines are 
included with the user's instructions to illustrate possible 
problem areas and the magnitude of the approximation associated 
with the model. 

Compressibility effects .- The selection of a compressibility 
correction scheme must take into consideration the configurations 
of interest and the calculation procedure. For example, vortex 
shedding and tracking is very dependent on cross sectional shape, 
and numerical transformations add significantly to computation 
time; therefore, it is advantageous that the compressibility 
transformation used have no effect on the cross sectional 
shape. With this guideline and based on similar experiences and 
reguirements (Ref. 12), a transformation which modifies only the 
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axial coordinate is needed. The Gothert Rule, described in Ref- 
erence 13 and others, is the choice for the compressibility 
transformation. A brief description of the method included in 
VTXCLD follows. 

The transformation from the compressible (x,y,z) coordinate 
system to the incompressible (x' $ y',z f ) system is 

x' = x/S 

y' = y 

z ' = z ( 18 ) 


where 


3 



(19) 


As a result of Equation (18), the body slope and angle of attack 
become 


dr' 

■ 3 

dr 

(20) 

dx ' 

dx 

a 1 

= 8a 


(21) 


and the velocity fields are related by 

u = u'/e 2 

v = v'/3 

w = w'/3 (22) 
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The vortex shedding scheme in VTXCLD requires that the compressi- 
bility correction be applied in a manner slightly different from 
the usual. At the initiation of the calculation, the actual 

missile body is transformed to the incompressible shape by the 

above stretching procedure, and the source distribution is ob- 

tained as described above. The modified incompressible flow 
conditions are used to calculate the velocity field associated 
with the body. This velocity field is transformed back to the 

compressible domain so that it is available for use in vortex 
tracking and compressible pressure coefficient calculations. The 
compressible pressure coefficient is used to locate separation 
and define the shed vorticity. Since the compressibility trans- 
formation has no direct effect on the cross sectional shape or 
the discrete wake vortices, the separation and vortex tracking 
calculations take place in the compressible domain. No compres- 
sibility effects are considered for the vortex induced velo- 
cities . 


Vortex Shedding Model 

The body vortex shedding model described in this section is 
nearly the same as the subsonic model presented in Reference 2; 
therefore, only the differences between the two models will be 
discussed in detail. 


Equations of motion .- The equations of motion of a shed nose 
vortex in the presence of other free vortices in the vicinity of 
a body in a uniform stream follow. In the circle (v) plane, the 
position of a vortex, r , is 


v 

m 


r 


m 


+ ix 


m 


and an image of r is located at 

m 


(23) 
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( 24 ) 


to satisfy the flow tangency condition on the body surface. 

In the real plane, the position of the vortex r is 

m 


0 = y m + iz 

m j, m in 


(25) 


The complex potential in the real plane is 


W( o) = $ — iip 


(26) 


and the corresponding velocity at r is 

m 


dW (a) . 

v - iw = — ^ [W (a)] ^1 
m m do dv m do , a=o 


(27) 


The complex potential of r is not included in Equation (27) to 

m 

avoid the singularity at that point. The derivative of the 
transformation is obtained from Equations (4) or (9). 


The total velocity at r in the crossflow plane is written 

m 


v - iw 

= G +G+G + G +0^ + 0 

V a g n m T r 

00 


(28) 


where each term in Equation (28) represents a specific velocity 
component in the o-plane. The first term represents the uniform 
flow due to angle of attack 


G = -i sina 1 + 
a 


dv | 

dal a= a 


(29) 
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The second term represents the uniform flow due to angle of yaw. 


Gg = -sing 


1 



d v | 

da ' o=a 

m 


(30) 


In compressible flow, the velocity components from Equations (29) 
and (30) used in the surface pressure calculations include ef- 
fects of the compressibility transformation described in a pre- 
vious section. 


The next term represents the influence of all vortices and 
their images, with the exception of r . 


N 


; n = 1 l 


n 


i 2 tt r V 
n = l o « 


m 




(v » /r o ,-(v n /r o ) 


The next term is due to the image of r . 

m 


dyi 

dal 
n*m a~o 


m 


(31) 


m 


G m 1 2irr V 
o 00 




dv 

da 


o = o 


m 


(32) 


The fifth term in Equation (28) represents the potential 


of T m in the o-plane and is written as 


G T ■ -* W- [|] h S3 (33) 

00 a = a^ 

m 

The last term in Equation (28) represents the velocity components 
induced by the body source singularities representing the volume 
effects. 
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v - iw 
G r = V 


( 34 ) 


where v r and w r are components of v R from Equation (11). 

Since the body source singularities are three-dimensional, 
they contribute an induced axial velocity, u r , given by Equation 
(10). Compressibility effects are included through the stretch- 
ing of the axial coordinate and subsequent transformation of the 
induced velocity components. 


The differential equations of motion for r are 

m 


da 


m 


v - iw 
m m 


dx 

V cos a + 

U 


OC Q 

r 

where 

0 = 

y - iz 


m 

m m 


Therefore , the two equations 

which must 

be 

body length to determine the 

trajectory of r 

dy 



dx 

_ 

V cos a + 

u 


oc C 

r 

and 

dz 

w 


m 

m 


dx 

V cos a + 

u 


oc c 

r 

There are a pair of equations like (37) 

and 


(35) 


(36) 


m 


(37) 


(38) 


in the field. As new vortices are shed, the total number of 
equations to solve increases by two for each added vortex. These 
differential equations are solved numerically using a method 
which automatically adjusts the step size to provide the speci- 
fied accuracy. 
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Surface pressure distribution .- The surface pressure dis- 
tribution on the body is required to calculate the forces on the 
body and the separation points. The surface pressure coefficient 
is determined from the Bernoulli equation in the form 



(39) 


(40) 


(41) 


where U is the total velocity (including V ) at a point on the 

oo 

body. 


The last term in Equation (41) represents the axial velocity 
missing from the two-dimensional singularities in the flow 
model. In this case, the shed vortices and the doublet part of 
Equations (29) and (30) are the only singularities contributing 
to this term. The complex potential of these two-dimensional 
singularities is 


W( v ) = <|> + i\|> 


n r 

= y _H 

2 u 

n+1 


+ l 


-i Jtn(v - 


V sina - 


v ) + i Hn 
n 


— V sing 

V or 


n 


(42) 
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The "unsteady" term in Equation (41), evaluated on the body 
surface, is 


= Real 221»i 
dx dx 


r=r 


(43) 


Equation (43) becomes 


, x nt 

g <|? _ y n 

dx ^ 2 it 

n = l 
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r dx 
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, X 2 
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(Xr 
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-5 r d X dx , dx dr 1 

X n r o ) L 2xX n dx - + 2lT n dx r o dx 2 r o T n dx~J 

2 2 
2 2 c o o 

(Tr n " T n r o } + Ur n ‘ Vo> 


2 V dr. 


r dx 
o 


( X sina - x sing) 


(44) 


where 


2 2 2 

r = x z + X z 

n n n 


(45) 


Separated wake »- The separated wake on the lee side of the 
body is made up of a large number of discrete vortices, each 
vortex originating from one of two possible locations at each 
time step in the calculation. The major portion of the lee side 
vortex wake has its origin at the primary separation points on 
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each side of the body. The remainder of the wake originates from 
the secondary separation points located in the reverse flow re- 
gion on the lee side of the body. Both of these points are il- 
lustrated in the sketch of a typical crossflow plane of an ellip- 
tic cross section body shown in Figure 3. The mechanics of the 
calculation of the individual vortices follows. 

As described in References 1 and 2, the pressure distribu- 
tion for the primary flow in the crossflow plane is referenced to 
the conditions at the minimum pressure point, and a virtual ori- 
gin for the beginning of the boundary layer is calculated. The 
adverse pressure distribution downstream of the minimum pressure 
point is considered with either Stratford's laminar (Ref. 14) or 
turbulent (Ref. 15) separation criterion to determine whether or 
not separation has occurred. These criteria, based on two- 
dimensional, incompressible, flat plate data, are adjusted for 
three-dimensional crossflow effects in Reference 2. The laminar 
separation criterion states that the laminar boundary layer sepa- 
rates when 




sina 

c 


In a turbulent boundary layer, separation occurs when 


(46) 


C 


P 



(Re^x 10 


- 0.1 


0.35 sina 

c 


(47) 


The sina c in Eguations (46) and (47) is the three-dimensional 
modification. 

If the criteria indicate a separation point, the vorticity 
flux across the boundary layer at separation is shed into a 
single point vortex whose strength is 
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I 


I 




1 assuming no slip 

described below. 

I 

The initial position of the shed vortex is determined such 
that the surface velocity in the crossflow plane at the separa- 
tion point is exactly canceled by the shed vortex and its 
image. When this criterion results in a vortex initial position 
that is too near to the body surface, certain numerical problems 
can cause difficulty in calculating the trajectory of this vor- 
tex. If the initial position of the vortex is nearer than five 
percent of the body radius from the surface, the vortex is gener- 
ally placed five percent of the equivalent body radius from the 
surface. This position represents the approximate thickness of 
the boundary layer. 

The calculation of secondary separation is carried out in 
the same manner as described above for primary separation. It is 

i 

i necessary that a reverse flow region exist on the lee side of the 

• body and that a second minimum pressure point be found in this 

region. Surface flow visualization of secondary separation re- 
gions on bodies at high angles of attack in subsonic flow are 
shown in Reference 16. For purposes of this analysis, the re- 
verse flow is assumed to be laminar from the lee side stagnation 
point to the secondary separation point, and Stratford's laminar 
criterion is used to locate secondary separation. Laminar sepa- 
ration in the reverse flow region is expected because of the low 
velocities on this portion of the body. The vortex released into 
the flow at the secondary separation point has the opposite sign 
of the primary vortex and is generally weaker in strength. 

Transition . - In general, the primary separation on a missile 
is either laminar or turbulent, depending on the Reynolds number 


u. 


_ _ Ax 

_„2 cos a 
« 2V c 


(48) 


at the wall. 6 is the vortex reduction factor 
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of the flow; however, as noted in Reference 17, there is an angle 
of attack and Reynolds number range over which the character of 
the separation is transitional. Flow visualization of the tran- 
sition phenomenon on missile shapes is available in References 16 
and 18. Body vortex shedding technology is not developed to the 
level that prediction of the beginning of transition nor the 
length of the transition is practical. Correlation of the lim- 
ited amount of data available has permitted the authors to define 
approximately the limits of transition. The results of the cor- 
relation indicate that the separation switches from laminar to 
turbulent when the Reynolds number based on the length from the 
nose times the sine of the angle of incidence is approximately 
7x10^. The length of the transition region is about one body 
diameter. More analysis is reguired to verify the correlation, 
and the authors recommend that transition effects be included 
with great care. 

Reduction factor .- The empirical vortex reduction factor, 

6, in Equation (48) is widely used in discrete vortex methods to 
provide better agreement between measured and predicted vorticity 
strength in the wake. Various investigators (e.g., Refs. 2, 4, 
and 19) recommend a value in the range 0.6 < 6 < 1.0 for subsonic 
flow with 0.6 being the more common value. Based on the experi- 
ence of the authors, the proper value of <5 can vary with the 
geometry and flow conditions. For example, 6 = 0.6 is an appro- 
priate value for an ogive-cylinder missile in subsonic flow, 
but 6 = 1.0 is a better value for a boattail body in subsonic 
flow. Unless otherwise noted, 6 = 0.6 is used for all predicted 
results presented in this report. 

Forces and moments .- The forces and moments on the body are 
computed by integration of the pressure distribution around the 
body. At a specified station on the body, the normal-force coef- 
ficient on a Ax length of the body is 
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where r' is the distance from the axis of the body to the body 
surface and e' is the local slope of the body in the crossflow 
plane. This is illustrated in the sketch in Figure 4. For cir- 
cular bodies, r' = r Q and 3 * = 0 - The total normal force coeffi- 
cient on the body is 


_ N Dr. 

C., = — — = — I c dx 
N q S S J n 

- oc o 


(50) 


and the pitching-moment coefficient is 
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The center of pressure of the normal force is 

x cp x ° 

n m 


*ref *ref 


(52) 


Similarly, the side-force coefficient on a Ax length of the 
body is 


r As i o 

I A - 1 2 IT 

= L * J = - i / C r' sing ' dg 
q D D J p 


The total side-force coefficient on the body is 


F D L 

C = ^ / c dx 

y q S S * y 


and the yawing-moment coefficient is 



The center of pressure of the side force is 


(53) 


(54) 
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A noncircular body at an arbitrary roll angle may experience 
a rolling moment caused by the nonsymmetry of the loading around 
the body. The total rolling moment on the body is calculated by 
summing the moments of the individual components of normal force 
and side force around each cross section and integrating over the 
body length. Using Eguations (4) and (5), the rolling moment 
coefficient is calculated as 
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Though the primary purpose of the vortex shedding method is 
not to predict drag or axial force on the missile, a procedure to 
estimate both pressure drag (excluding base drag) and skin fric- 
tion is included. The pressure contribution to the axial-force 
coefficient is 



, L 2u 

- 5 J / C, 


dr' 

dx 


dp dx 


(58) 


The local skin-friction coefficient, based on the assumption of a 
l/7th power law velocity profile in the boundary layer, is 


C 


f 


.0592 



(59) 


which produces a drag coefficient due to friction 


’ = — — f r C dx 
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(60) 
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At high angles of incidence, streamlines around the body are 

inclined at approximately 

a = tan* (2 tan a ) (61) 

s c 

therefore, the axial component of the friction drag is 

C = C n cos a (62) 

A f D f s 

The total axial-force coefficient is the sum of Eguations (58) 
and ( 62 ) . 

Vortex core . - The diffusion core model (Ref; 2) for the 
point-vortex induced velocities removes the singularity at the 
vortex origin and effectively reduces the velocities near the 
vortex. The tangential velocity induced by a single point vortex 
is written as 



where r is the distance from the vortex to the field point and x 
is a measure of the age of the vortex. The induced velocity from 
this core model is illustrated in the following sketch. 



Sketch - Vortex Induced Velocity 
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The vortex core model represented by Equation (63) has re- 
ceived considerable attention in the context of body vortex in- 
duced effects, and it has a number of shortcomings. Since the 
exponential term is a function of r, the flow medium (v), and the 
age of the vortex (x), the core is constantly changing size as 
the vortex moves through the field. Under certain conditions, 
the core radius, denoted as r c in the sketch, can become very 
small and the induced velocity becomes unrealistically large. In 
an attempt to further modify the core model to keep the induced 
velocities within physically realistic limits, the following 
modification was made in Reference 3. The location of the maxi- 
mum induced velocity, r c in the sketch, is fixed at a specific 
radius to be selected by the user. Given a core radius, the 
vortex induced velocity is 


1 - exp 

The core model described by Equation (64) is included in the 
version of the code described in this report, and in a later 
section, guidelines for the selection of an appropriate core 
radius are presented. Results indicate that this simple core 
model provides adequate smoothinq for the necessary velocity 
calculations . 



-[ 


1.256 


r^ 

r 2 J 
c 


(64) 


Other investigations of discrete vortex models must use a 
core model of some type. Although there are many different core 
models, all serve the same purpose as those described above, and 
nearly all are directed at eliminating the singularity and reduc- 
ing the maximum induced velocity. The code VTXCLD described in 
this report is easily modified to incorporate another core model 
if the user desires to do so. 

Asymmetry .- An axisymmetric missile at large angles of at- 
tack in subsonic flow will exhibit asymmetries in the vortex 
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field and associated induced side force. The vortex cloud ap- 
proach described herein will predict no asymmetry at any angle of 
attack on a missile with a vertical plane of symmetry unless the 
solution is perturbed. Every investigator (e.g. f Refs. 2, 5, and 
19) has a favorite means of introducing a perturbation to the 
solution, and no good argument can be made for one over an- 
other. Unfortunately, the size and type of perturbation has an 
effect on the final solution; therefore, some empircal informa- 
tion is reguired to calibrate the perturbation for a particular 
geometry and flow condition. 

In Reference 2 and in VTXCLD, the perturbation is introduced 
by a small modification of the predicted symmetric separation 
points on either side of the missile at the first few integration 
steps near the nose. The solution is then allowed to develop on 
its own. Some results for various magnitude perturbations are 
presented in Reference 2 and this report. 

PROGRAM VTXCLD 

Program VTXCLD is a rational flow model which predicts the 
nonlinear forces and moments acting on slender circular and non- 
circular bodies at high angles of attack in steady subsonic flow 
conditions. The nonlinear effects are obtained by modeling the 
lee side separation vorticity with discrete vortices in the 
crossflow planes. The code has application as an engineering 
prediction method for use in preliminary design and analysis 
studies of missiles at high angles of incidence, and as such, it 
is important that the code be simple and economical to use and 
reguire a minimum of empirical information. However, it is 
egually important that the code adeguately represent the impor- 
tant nonlinear flow phenomena and simulate the major physical 
features of the flow field. 
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The following sections are a user's manual for program 
VTXCLD. These sections include descriptions of the program, the 
calculation procedure, input preparation, sample input and output 
cases which illustrate various program features, and program 
limitations . 


General Description 

This section provides a general description of program 
VTXCLD and its various subroutines. The code is written in modu- 
lar form so that as improvements in flow models become available 
or as other modifications are required they are easy to incorpo- 
rate. The flow of the calculation is described in this report 
and detailed comments are provided throughout the code to assist 
the user in understanding the order of calculation. 

The computer code consists of the program VTXCLD and 39 
subroutines. The overall flow map of the program is shown in 
Figure 5 where the general relationship between the subroutines 
and external references can be seen. Communication between the 
program modules is handled primarily by named common blocks. A 
cross reference table for the calling relationship between the 
program subroutines is shown in Figure 6. A similar table for 
the named common blocks is shown in Figure 7. 

The program is written in standard FORTRAN. Execution on a 
VAX 11/750 can vary from 10 seconds to more than 45 minutes de- 
pending on various factors such as the geometry, the integration 
interval (DX), the number of shed vortices, the use of flow sym- 
metry, and the flow conditions. 
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Subroutine Description 


I 

I 

This section briefly describes the main program VTXCLD and 
I its various subroutines. 

I 


VTXCLD This is the main program for calculating the forces 
and moments on a slender body in a steady flow 
condition including the nonlinear effects of lee 
side separation vorticity. 

Transforms a specified point in the circle plane to 
its corresponding point in the body plane. 

Calculates velocity term due to transformation for 
noncircular bodies. 

Sets up a table of points on the actual body cor- 
responding to specified points on the circle. 

Organizes the flow through VTXCLD. It organizes 
calls to predict vortex shedding, to calculate 
surface velocities and pressures, and to calculate 
the overall forces and moments on the body. 

Transforms a specified point in the body plane to 
its corresponding point in the circle plane. 

COMBIN Combines vortices which are separated by a distance 
less than or egual to " RGAM " . (The resulting vor- 

tex strength is the sum of the individual strengths 
and its location is the centroid of the combined 
vortices ) . 
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CONFOR 


CP POT 


DFEQKM 


DPHIDT 


DZDNU 


F 


FDPSDR 


FORCE P 


FPSI 


FPVEL 


Calculates the coefficients of the numerical con- 
formal mapping function. Transforms the region 
outside of a polygonal shape with a vertical plane 
of symmetry to the region outside of a circle. 

Calculates a potential pressure distribution at a 
given x station (the effects of separation vortices 
are not present in the C p calculation). This is 
used in determining a turbulent separation point 
for the transition region. 

Obtains vortex paths between stations x and x to Dx 
using a Kutta-Merson integration scheme. 

Calculates the two-dimensional unsteady pressure 
term for use in the unsteady Bernoulli eguation. 

Calculates the differential of the transform of a 
polygonal shape to a circle. 

Calculates the derivatives of the vortex eguations 
of motion which are used in conjunction with DFEQKM 
to obtain the vortex paths. 

Calculates the axial velocity due to a three- 
dimensional source distribution. 

Calculates the forces and moments on the body by an 
integration of the surface pressure distribution. 

Calculates the body shape represented by a three- 
dimensional source distribution. 

Calculates the velocity at specified field points. 
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INPT 

INVERS 

PLOT 

RDRDX 

RDRDXI 

SEPRAT 

SHAPE 

SKIN 

SORDIS 

SORVEL 


Reads in and prints out the body geometry, the flow 
conditions, various program options, and restart 
guantities . 

Solves simultaneous equations. Ax = B. 

Creates a line printer plot of the body and its 
vortex wake at any given x station. It is composed 
of the following subroutines: PLOTV, PL0TA2 , 
P LOTA 5 , PL0TA6, PL0TA7 , and PL0TA8 . 


Calculates r and dr/dx at a given x station. 

Initializes the calculation of r and dr/dx for 
subroutine RDRDX. 

Predicts the position in the crossflow plane where 
separation occurs based on the Stratford criteria 
for laminar and turbulent flows. 

Determines the body geometric parameters from a 
table look-up of input data. 

Estimates the axial skin friction force on the 
body . 

Calculates the source distribution representing an 
axisymmetric body of revolution given the body 
shape and slope at each x station. 

Calculates perturbation velocities due to a three- 
dimensional source distribution. Used by subrou- 
tine SORDIS in calculation of the source distribu- 
tion . 
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SRCVEL 


Calculates perturbation velocities due to a three- 
dimensional source distribution. 


SUM 

TRANSN 


VCENTR 


Analytical mapping of an elliptical cross section. 

Calculates the body separation points in the tran- 
sition region. 

Calculates the centroid of each vortex field. 


VLOCTY Calculates the velocity components in the cross 

flow plane at a specified field point. 

VTABLE Sets up the arrays of free vorticity for two 

cases: (1) The vorticity in the working array GAM 

is distributed into the individual arrays GAMP, 
GAMM, GAMR , GAMA, or (2) the vorticity in the indi- 
vidual arrays GAMP, GAMM, GAMR, GAMA is distributed 
into the working array, GAM, in the proper order. 


Z Calculates the transformation of a polygonal shape 

to a circle of radius RZ . 

Program Limitations and Suggestions 


Program VTXCLD predicts the aerodynamic forces acting on 
slender bodies alone at high angles of attack, with some limita- 
tions. Certain of the limitations that follow are basically 
rules-of-thumb for the use of the program which have been deter- 
mined through program verification and experience with the code. 


Mach number . - The major limitation in the program is the 
Mach number or compressibility correction. The compressibility 
correction added to VTXCLD is based on small perturbation theory 
which restricts the body to be slender, the angle of attack to be 
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small , and the Mach number to be approximately 0.8 or less. The 
slender body restriction is not a new restriction as this assump- 
tion is inherent in the crossflow analogy used by VTXCLD to pre- 
dict lee side vorticity. These restrictions result from the 
assumptions made for small perturbation theory, and any cases 
outside these restrictions should be viewed with caution; how- 
ever, VTXCLD has been applied to blunt nose shapes such as spher- 
oids or ellipsoids with good success. Typical missile shapes 
with slender noses pose no problems to the method. 

Source distribution .- The following suggestions should be 
observed when using the source distribution calculation in 
VTXCLD. The source distribution in VTXCLD places three-dimen- 
sional point sources on the axis of the body (or eguivalent body 
in a noncircular case). For subsonic Mach numbers it is neces- 
sary to close the body to calculate a smooth source distribu- 
tion. Generally, the body can be closed at its base by including 
a smooth boattail, and a successful technigue for most cases is 
to close the model with its nose. The source distribution does 
not handle rapid changes or discontinuities in body slope well. 

The source distribution calculation often has a problem with 
short pointed noses as well. Short blunted noses seem to cause 
no difficulty because they resemble Rankine bodies, a single 
point source in a uniform stream. Sharp pointed noses often 
cause an oscillation in the source strength on the afterbody 
resulting in an oscillation in the velocities and thus an oscil- 
lation in the axial pressure distribution. This has been ob- 
served with a two-caliber ogive nose but was not apparent for a 
2.5-caliber ogive. The oscillations were minor in their effect 
on the total calculation and result, but the user should be aware 
of their possibility. 
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Incidence angle .- The well-known and little-understood phe- 
nomenon of wake asymmetry at large incidence angles (a > 30°) is 

c 

available as an option in VTXCLD. As described in a previous 
section, the naturally symmetric nature of the vortex shedding 
analysis must be perturbed to produce an asymmetric solution; 
however, the asymmetric capability within VTXCLD has been demon- 
strated and will be described in the RESULTS section. The user 
must be aware of the limitations associated with this portion of 
the code in order to apply appropriate judgement in the evalua- 
tions of the results. Based on the author's experience with 
asymmetry, care must be taken in the choice of the perturbation 
parameter and its effect on the results. This is to serve as a 
warning, and the user should be aware that this portion of the 
code is in need of additional development. 

Transition .- The transition capability included in the code 
is also in an early stage of development. The character of the 
vortex wake, whether originating from laminar or turbulent sepa- 
ration points, has an important effect on the predicted nonlinear 
aerodynamics characteristics. As previously discussed, addi- 
tional development is needed to better understand the origin and 
length of the transition regions. The user is cautioned to be 
very careful in the interpretation of predicted results in flow 
regimes which may be in a transition region as defined by Lamont 
in Reference 17. 

Secondary separation .- The capability to predict secondary 
separation or separation in the reverse flow region on the lee 
side of the body is included in VTXCLD. As discussed in a later 
section, there are certain flow conditions in which secondary 
separation has an important role; for example, asymmetric shed- 
ding at large incidence angles. Numerical problems associated 
with tracking the secondary vorticity have been observed, and 
some additional effort is reguired to make secondary separation a 
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useful option. Therefore, it is recommended that this option be 
■ treated as a research feature and not used for production calcu- 
lations. 

Error Messages and Stops 

j The code VTXCLD has numerous internal error messages. These 

messages are generally explicit and are described below. There 
are several execution terminations at numbered stops within the 
program. These are described in the table below with some sug- 
J gestions to correct the problem. 

The error message "CONVERGENCE NOT OBTAINED IN INTEGRATION" 
j is printed out when the trajectory calculation does not converge 

within the specified tolerance, E5 . Also printed out is the 
axial station, XI, and the integration interval, H. This error 
i results in a STOP 40 which is discussed in the table to follow. 

I 

! 

In the calculation of the source distribution it is possible 
to get the message "MATRIX IS SINGULAR. " This occurs when one of 
the pivot elements in the influence matrix is zero. This error 
results in a STOP 300. 

The other major error messages are found in the numerical 
conformal mapping routine CONFOR. The message "***ERR0R IN THE 
SUM OF EXTERIOR ANGLES***" is printed if the routine has diffi- 
culty calculating appropriate body angles within a specified 

tolerance. Another message "***COEFFICIENT AN(1) = GREATER 

THAN +- 0.0001 AFTER ITERATIONS OUT OF " is printed on 

the rare occasions when the mapping procedure is having diffi- 
culty converging on a solution. These messages are informational 
and do not result in a program STOP. However, the user is ad- 

vised to check the transformed shape carefully with the actual 
shape before using the results. 
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Termination 

Message 

STOP 1 
STOP 40 


STOP 41 
STOP 44 


Description 

This is the normal stopping condition. Pro- 
gram VTXCLD 

The vortex trajectory integration has failed 
to converge at an axial station. This problem 
is usually caused by: (1) two or more vor- 

tices in close proximity rotating about one 
another so rapidly that the integration sub- 
routine cannot converge on a solution, or (2) 
a vortex is too near the body surface. These 
situations are often dependent on specific 
flow conditions and the body shape. The fol- 
lowing steps should be tried to resolve the 
problem. increase the integration tolerance 
E5, or increase the size of the combination 
parameter RGAM so that it is greater than the 
distance between the troublesome vortices. It 
often helps to increase the integration inter- 
val DX and restart the calculation procedure 
from a previous station. Subroutine BODY 

The vortex trajectory calculation has resulted 
in a vortex inside the body. This is a rare 
but serious problem which indicates a major 
problem with the solution. Subroutine BODY 

The total number of shed vortices exceeds 
200. To correct this combine some vortices 
and restart the calculation, or increase the 
integration interval, DX, so that fewer vor- 
tices are shed. Subroutine BODY 
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STOP 45 


The maximum number of positive or negative 
primary separation vortices exceeds 70. Use 
the same remedy as for STOP 44. Subroutine 
BODY 

; 

STOP 46 The maximum number of positive or negative 

I secondary (reverse flow) separation vortices 

exceeds 30 . Use the same remedy as for STOP 
44 but only combine like vortices. If one of 
the above steps does not work, the reverse 
flow calculation can be turned off by setting 
NSEPR = 0. Subroutine BODY 

STOP 201 There is an error in the vortex tables in 

STOP 202 subroutine VTABLE . Check the input values of 

NBLSEP , NSEPR, NVP, NVM , NVA, NVR, and RGAM in 
Items 1 through 6. 

STOP 300 There is a singular influence matrix in the 

solution for the three-dimensional source 
distribution. Check the input geometry for 

errors. Subroutine INVERS 

STOP 301 The maximum number of iterations is exceeded 

in the source calculation. Check the input 
geometry. This problem can occur if the body 
has a long slender nose followed by a cylin- 
drical section and then a long slender tail. 
If this is the case, one solution is to close 
the body off with a blunter and/or shorter 
tail. Subroutine SORDIS 
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Input Description 


This section describes the input data and format required by 
program VTXCLD. The input formats are shown in Figure 8, and the 
definitions of the individual variables are provided in the fol- 
lowing paragraphs. Note that some length parameters in the input 
list are dimensional variables; therefore, special care must be 
taken that all lengths and areas are input in a consistent set of 
units. 


The remainder of this section includes a description of each 
of the input variables and indices reguired for the use of pro- 
gram VTXCLD. The order follows that shown in Figure 8 where the 
input format for each item and the location of each variable on 
each card is presented. Data input in I-format are right justi- 
fied in the fields, and data input in F-format may be placed 
anywhere in the field and must include a decimal point. 


Item 1 is a single card containing 15 integer flags. These 
flags are right justified in a five column field and are used to 
specify several options available in VTXCLD. 


Item Variable Description 

1 (1515) NCIR Cross section flag, determines the 

type of body cross section. 

= 0, circular body 
= 1, elliptical body 
= 2, arbitrary body 


NCF - Mapping flag, specifies whether the 

reguired mapping is input or calcu- 
lated. (NCIR = 2) 

= 0, calculate mapping 
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= 2, input mapping 


ISYM 


NBLSEP 


NSEPR 


NDFUS 


NDPHI 


INP 


Vortex shedding symmetry flag. 

= 0, symmetric shedding 
= 1, nonsymmetric shedding 

Separation flag. 

= 0, no separation 
= 1, laminar separation 
= 2 , turbulent separation 

Reverse flow separation flag. 

= 0, no reverse flow separation 
(preferred) 

= 1, reverse flow separation 

Vortex core model flag. 

= 0, potential vortex model 
= 1, diffuse vortex model (pre- 
ferred) 

Two-dimensional unsteady pressure 
term flag. 

= 0, omit 3$/3t term in Cp equation 
= 1, include 3<f>/3t term (preferred) 

Nose force flag. 

= 0 f estimate nose forces ahead of 
the first axial station (pre- 
ferred) 

= 1, zero nose forces ahead of the 
first axial station 
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NXFV Number of x stations at which the 

flow field is calculated and/or spe- 
cial output is generated. 

(0 < NXFV < 8) 

NFV Number of field points for flow 

field calculation. 

(0 < NFV < 200) 

NVP Number of +r vortices on +y side of 

body for a restart calculation. 

(0 < NVP < 70) 

NVR Number of reverse flow vortices on 

+y side of body for a restart calcu- 
lation. (0 < NVR < 30) 

NVM Number of -r vortices on -y side of 

body for a restart calculation. 

(0 < NVM < 70) 

NVA Number of reverse flow vortices on 

-y side of body for a restart calcu- 
lation. (0 < NVA < 30) 

NASYM Asymmetric vortex shedding flag. 

= 0 , no asymmetry 
= 1 , forced asymmetry 

Item 2 is a single card containing 10 integer flags/parame- 
ters. These values are used to specify several options available 
in VTXCLD . 
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Item 


Variable 


Description 


2 (1015) 


NHEAD Number of title cards in Item 3. 

(NHEAD > 0) 

NPRNTP Pressure distribution print flag. 

= 0, pressure distribution output 

only at special output sta- 

tions ( NXFV > 0) 

= 1, pressure distribution output 

at each x station 


NPRNTS Vortex separation print flag. 

=0, no separation output 

= 1, separation point summary (pre 
ferred) 

= 2, detailed separation calcula- 
tion (for diagnostic purposes 
only) 

NPRNTV Vortex field summary output flag. 

= 0, vortex field is output at 

special output stations only 
(NXFV > 0) 

= 1, vortex field is output at 
every x station 

NPLOTV Vortex field printer-plot flag. 

= 0, no plots 

= 1, plot full cross section, con- 
stant scale 

= 2, plot half cross section, con- 
stant scale 
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NPLOTA 


NTH 


NCORE 


NSKIN 


NCOMP 


= 3, plot full cross section, vari- 
able scale 

Additional output flag. 

= 0, no additional output 
= 1, additional output only at 

special x station (NXFV > 0) 

= 2, additional output at each 
station 

Numerical conformal mapping angle 
flag . 

= 0, calculate circle 0's for map- 
ping at 5° increments (pre- 
ferred ) 

= 1, input Q's 

Vortex core radius flag. 

= 0, maximum allowable local vortex 
core size is .05*SD 
= 1, no upper limit on the local 

vortex core size, the vortex 
core size is given by: RCORE*D 

Axial skin friction estimation flag. 
= 0, no skin friction estimation 
= 1, axial skin friction estimated 

Compressibility correction flag. 

= 0, incompressible calculation 
= 1, compressibility correction 

applied (Mach number must be 
input in Item 6) 
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Item 3 is a set of NHEAD title/summary cards which are 
printed at the top of the first output page. They are reproduced 
in the output just as they are input. 

Item Variable Description 

3 (20A4) TITLES NHEAD title/summary cards to be 

output on first page. 

Item 4 is a single card containing reference information 
needed to form the aerodynamic coefficients. 


Item Variable 

4 (5F10.5) REFS 

REFL 

XM 

SL 

SD 

Item 5 is a single card 
Item Variable 

5 (3F10.5) ALPHAC 

PHI 


Description 

Reference area, S re £. REFS >0. 

0 

Reference length, ref. REFL > 0. 

Axial position of the moment center, 
measured from the body nose. 

Body length, L. 

Maximum body diameter, D. 

containing the flow conditions. 

Description 

Angle of incidence (degrees). 

(0° < a < 90°) 
c 

Angle of roll, <J>, (degrees). 
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RE 


Reynolds number based on maximum 
body diameter, SD, and V .(V D/v) 

oc oc 


Omit Item 6 if NCOMP = 0. 

Item 6 is a single card which contains the freestream Mach 
number M . This item is included only if NCOMP in Item 2 is 

oc 

nonzero . 

Item Variable Description 

6 (F10.5) MACH Freestream Mach number. 

Item 7 is a single card which specifies the axial extent of 
the run, the transition region, and parameters associated with 
the vortex wake. 

Item Variable Description 


7 (8F10.5) XI 


XF 


DX 


XTRl 


Initial x station, XI > 0, dimen- 

sional quantity. (Note: XI > DX/2) 

Final x station, XF > XI, dimen- 
sional quantity. 

x increment for vortex shedding 
calculation. Typical value, 

DX = D/2. 

Beginning of transition region, 
dimensional quantity. 


XTR2 End of transition region, dimen- 

sional quantity. 
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EMKF 


Minimum distance of shed vortex 
starting position from the body 
surface . 

= 1.0, shed vortices located ac- 
cording to stagnation point 
criterion 

> 1.0, minimum radii position of 
shed vortices. Typical 
value, EMKF = 1.05 


RGAM Vortex combination factor. 

= 0.0, vortices are not combined 
( preferred ) 

> 0.0, radial distance within which 
vortices are combined 


VRF Vortex reduction factor, 6^ to 

account for observed decreases in 
vortex strength. 

= 0.6, for bodies with bases 
= 1.0, for closed bodies 


Item 8 is a single card 
ance, asymmetry perturbation 
criteria, and the core radius. 


containing the integration toler- 
information, alternate separation 


I tern Variable Description 

8 (8F10.5) E5 Error tolerance for the vortex tra- 

jectory calculation. Typical range, 
E5 = 0.01 to 0.05 or larger if ap- 
propriate . 
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XTABL 


Location at which a table of points 
on the body and on the circle are 
printed . 

= 0.0, no table of points is 
printed (preferred) 

> 0.0, table of points for all 
x < XTABL 

XASYMI Initial x station for asymmetric 

perturbation . 

XASYMF Final x station for asymmetric per- 

turbation . 

DBETA Angular displacement of separation 

points for asymmetry perturbation. 

SEPL Stratford laminar criterion. 



li 

o 

• 

o 

> 

program uses 
( preferred) 

SEPL = .087 


V 

o 

• 

o 

program uses 

input value 

SEPT 

Stratford turbulent 

criterion. 


11 

o 

• 

o 

■> 

program uses 
( preferred ) 

SEPT = .350 


o 

• 

o 

A 

program uses 

input value 

RCORE 

Ratio 

of the local 

core radius 


the local body diameter. Default 
value is .025. Maximum allowable 
value is .025*D/r Q . 

Item 9 is a single card containing two flags for the source 
distribution. The source distribution may be input (NSOR > 0), 
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calculated from the body geometry (NSOR = -1 ) , or omitted from 
the solution (NSOR = 0) . If the Mach number is not equal to zero 
(NCOMP > 0) then a different source distribution is required for 
each Mach number. 

Item Variable Description 

9 (215) NSOR Source distribution flag. 

= -1, calculate the source distri- 
bution from the input geome- 
try input 

= 0, no sources (not recommended) 

> 0, number of sources to be 

input 

NPRT Source distribution print flag. 

= 0, output source locations and 

strengths 

= 1, above output plus input and 

source geometry comparison 

= 2, above output plus source- 

induced surface velocities 

Omit Items 10 and 11 if NSOR < 0. 

Items 10 and 11 contain the nondimensional source locations 
and nondimensional source strengths, respectively. 

Item Variable Description 

10 (6E12.5) XSRC Nondimensional source locations, 

x/L. (NSOR values, 6 per card) 
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Item 


Variable 


Description 


11 (6E12.5) QS Nondimensional source strengths, 

o 

Q/4ttL V cos a . 

' op c 

Items 12 through 15 contain the body geometry information. 
This information must be included for all bodies whether in cir- 
cular or noncircular. 


Item Variable Description 

12 (15) NXR Number of entries in the body geome- 

try table. (1 < NXR < 101) 


Item Variable Description 

13 (8F10.5) XR Nondimensional axial stations in the 

geometry table, x/L. (NXR values, 8 
per card) 


Item Variable Description 

14 (8F10.5) R Nondimensional body radius at the 

NXR x/L-stations, r Q /L (NXR values, 
8 per card) . For an ellipse, r Q = 
(a + b)/2. For arbitrary bodies, r Q 
is equivalent radius. 


- 52 - 




Item 


Variable 


Description 


15 (8F10.5) DR Body slope at the NXR x/L stations 

dr Q /dx. (NXR values, 8 per card) 

Omit Items 16 and 17 if NCIR 1. 

Items 16 and 17 contain the nondimensional horizontal and 
vertical half axis lengths if the cross sections are elliptical, 
NCIR = 1. 


Item 


Variable 


Description 


16 (8F10.5) AE Nondimensional horizontal half axis 

length of the elliptical cross sec- 
tion at the NXR x/L-stations, a/L. 
(NXR values, 8 per card) 


Item Variable Description 


17 (8F10.5) 

BE 

Nondimensional vertical 

half 

axis 



length of the elliptical 

cross 

sec- 



tion at the NXR x/L-stations, 

b/L. 



(NXR values, 8 per card) 




Omit Items 18 through 29 if NCIR /= 2. 


Items 18 through 29 are used to specify a body of arbitrary 
cross section shape. This body should also be polygonal in shape 
(i.e., negative interior angle can cause problems in the numer- 
ical conformal mapping procedure). 
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Omit Item 18 if NTH = 0 


I tern 

18 (8F10. 


Item 

19 (215) 


Item 

20 (8F10 


Omit Items 


Variable Description 

5) TH Special table of circle angles for 

arbitrary body mapping, (i = 1, 73), 
or to be used only if the density of 
points on the cross section needs 
adjusting to improve the resolution 
of pressure distribution. Not rec- 
ommended under normal conditions. 


Variable Description 

MNFC Number of mapping coefficients for 

the arbitrary body. (MNFC < 100) 

MXFC Number of axial stations at which 

the arbitrary body is specified. 
(1 < MXFC < 10). For a similar 
shape at all cross sections, MXFC=1. 


Variable Description 

.5) XFC x stations at which the arbitrary 

body is specified. (MXFC values, 8 
per card). For a similar cross 
section body, MXFC=1, XFC(l) < XI. 

21 through 23 if NCF > 0. 
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Items 21 through 23 are repeated for each of the MXFC axial 
stations when the transformation coefficients are to be calcu- 
lated (NFC = 0) . 


Item Variable Description 

21 (15) NR Number of points describing the 

arbitrary body cross section 
(2 < NR < 30). (Note: NR must be 

the same for all x stations. ) 


I tern Variable Description 

22 (6E12.5) XRC Horizontal coordinates of the arbi- 
trary body cross section. (NR 

values, 6 per card). The convention 
for ordering the coordinates from 
bottom to top in a counter-clockwise 
fashion is observed. Figure 9. 
Right/left body symmetry is re- 
guired. 


23 (6E12.6) YRC Vertical coordinates of the arbi- 
trary body cross section. (NR 

values, 6 per card). See Figure 9. 


Omit Items 24 through 29 if NCF = 0. 


Items 24 through 29 are repeated for each of the MXFC axial 
stations if the numerical transformations are available from a 
previous calculation. The only purpose for providing the option 
to input the mapping parameters is to eliminate the need to re- 
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compute the numerical mapping in subsequent runs and save compu- 
tation time. 


Item 

24 (2F10.5) 


Item 
25 (15) 


Item 

26 (6E12.5) 


Item 

27 (6E12.5) 


Variable 

ZZC 


Description 

Mapping offset distance. 


RZC Equivalent radius of transformed 

cross section. 

Variable Description 

NR Number of points describing the 

arbitrary body cross section 

(NR < 30). (Note: NR must be the 

same for all x stations.) 


Variable Description 

XRC Horizontal coordinates of the arbi- 
trary body cross section. (NR 

values, 6 per card). See Figure 9 
and convention for Item 22 . 


Variable Description 

YRC Vertical coordinates of the arbi- 

trary cross section. (NR values, 6 
per card). See Figure 9. 
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I 

i 

I 

I 

Item Variable 


I 28 (6E12.5) THC 

I 

l 


Item Variable Description 

29 (6E12.5) AFC Mapping coefficients. (MNFC values, 

6 per card) 

Omit Items 30 and 31 if NXFV = 0. 

Item 30 is a single card containing the x stations at which 
additional output is printed and/or calculated if NXFV > 0. 

Item Variable Description 

30 (8F10.5) XFV x stations at which additional out- 

put is printed or calculated. (NXFV 
values, 8 values maximum) 

Omit Item 31 if NFV = 0. 

Item 31 is a set of NFV cards which contains the coordinates 
of the field points where the velocity components are to be cal- 
culated at each of the axial stations in Item 30. 

Item Variable Description 

31 (2F10.5) YFV,ZFV y and z coordinates of the field 

points for the velocity field calcu- 
lation. (NFV cards, maximum of 200) 


Description 

Circle angles 0, of body points. 
(NR values, 6 per card) 
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It is important that all points lie 
outside the body surface at all 
axial stations. 

Omit Items 32 through 37 if NVP + NVR + NVM + NVA = 0 

Items 32 through 36 provide information needed to restart 
the calculation procedure. It is not necessary to input all 
four types of vorticity in a restart calculation. 

Item 32 is a single card containing the nose force and mo- 
ment coefficients at the restart point; that is, X = XI . These 
values may be set equal to zero, but the resulting forces and 
moments from the current calculation will be in error. 


Item Variable 


32 (5F10.5) CN 

CY 

CA 

CDI 


CM 

CR 

CSL 


Description 

Normal-force coefficient 

Side-force coefficient 

Axial pressure force coefficient 

Axial skin friction force coeffi- 
cient 

Pitching-moment coefficient 
Yawing-moment coefficient 
Rolling -moment coefficient 
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Item 33 is a set of NVP cards which specify a cloud of posi- 
tive primary body separation vorticity, usually on the right side 
of the body. The variable XSHEDP associated with each vortex is 
used only to identify individual vortices and permit the user to 
follow individual vortices during the calculation. 


Omit Items 33 if NVP = 0. 


I tern 


Variable 


33 (4F10.5) GAMP 

YGP 

ZGP 

XSHEDP 


Description 

+y primary vorticity, r/V . 
Horizontal coordinate, y, of vortex. 
Vertical coordinate, z, of vortex. 
Origin of the vortex. 


Item 34 is a set of NVR cards which specify a cloud of re- 
verse flow (secondary) body separation vorticity. This array is 
a convenient way to put an arbitrary cloud of additional vorti- 
city which is to be maintained separately from the other body 
vorticity in the field. 


Omit Item 34 if NVR = 0. 


Item 

Variable 

Description 


34 (4F10.5) 

GAMR 

+y reverse flow vorticity. 

r/v . 
00 


YGR 

Horizontal coordinate, y, 

of vortex. 


nr»n 
£j \JJ\ 

Vertical coordinate, z, of 

vortex . 
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XSHEDR 


Origin of the vortex 


Item 35 is a set of NVM cards which specify a cloud of nega- 
tive primary body separation vorticity, usually on the left side 
of the body. In the case of a symmetric flow field (ISYM=0), 
this block of vorticity is automatically defined as the mirror 
image of the positive vorticity input in Item 33. 

Omit Item 35 if ISYM = 0 or NVM = 0. 


I tern 

Variable 

Description 

35 (4F10.5) 

GAMM 

-y primary vorticity, r/V . 

OP 


YGM 

Horizontal coordinate, y, of vortex. 


ZGR 

Vertical coordinate, z, of vortex. 


XSHEDR 

Origin of the vortex. 


Item 36 is a set of NVA cards which specify a cloud of re- 
verse flow (secondary) body separation vorticity. This block of 
vorticity is analogous to Item 34, and it is omitted if NVA = 0 
or if ISYM = 0. As with the previous item, this block of vorti- 
city is automatically defined as the mirror image of the negative 
vorticity input as Item 34 if symmetry is reguired by ISYM = 0. 

Omit Item 36 if ISYM = 0 or NVA = 0. 


I tem 

Variable 

Description 


36 (4F10.5) 

GAMA 

-y reverse flow vorticity. 

r/v . 

OP 


YGA 

Horizontal coordinate, y. 

of vortex. 
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ZGA 


Vertical coordinate, z, of vortex. 


ZSHEDA Origin of the vortex. 

Item 37 contains the pressure distribution at the initial x 
station for a restart calculation. If a symmetric case is being 
considered, ISYM=0, only half of the pressure distribution is 
needed. 


Item 


Variable 


37 (3F10.5) YBO 

ZBO 

CPO 


Description 

y body coordinate. 

z body coordinate. 

pressure coefficients at 
(YBO, ZBO). 


This concludes the input deck for a single run of VTXCLD. 
The code is not set up to stack input for multiple cases because 
of the sometimes long and generally unpredictable run times. 

Input Preparation 

In this section, some additional information is provided to 
assist the user in the preparation of input for certain selected 
problems. The previous section on the input description must be 
used to understand the individual variables which go into VTXCLD, 
and this section will permit the user to select the appropriate 
input to get optimum results from the code. 


Numerical mapping .- The specification of the appropriate 
numerical mapping parameters (Items 18 through 29) depends on the 
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tion time, the fewest possible terms in the transformation series 
( MNFC in Item 19) should be used which will produce an adequate 
mapping. Each different shape should be checked by the user to 
determine the proper number of coefficients. For example, a 
simple cross section which is similar to a circle may require as 
few as 10 coefficients. Others may require as many as 20 coeffi- 
cients for an excellent mapping. If a large number of calcula- 
tions are to be made for a noncircular shape, it is worthwhile 
for the user to try different numbers of coefficients in an ef- 
fort to find an optimum number. Not only should the shapes be 
compared carefully, the unseparated flow pressure coefficients 
should be compared to evaluate the quality of the mapping. 

The numerical mapping now included in VTXCLD has been used 
for a variety of cross section shapes, and very few problems have 
been observed. One general problem area is bodies with concave 
sides. The mapping procedure has difficulty converging on a set 
of coefficients. Usually the body shape can be modified 
slightly, to give the concave region a flat or slightly convex 
shape, and the mapping converges easily. This small change in 
the shape probably has very little effect on the predicted pres- 
sure distribution. This is an example of the type of problem to 
which the user should be alert. 

Integration interval .- Choice of the appropriate integration 
interval ( DX in Item 7) depends to a large extent on the type of 
results of interest. If the user is interested in detailed nose 
loads and separation on the nose, where the radius is changing 
rapidly, the internal should be smaller than the recommended 
D/2. On very long bodies, the recommended interval is usually 
adequate for accurate predictions. 

Vortex core .- The core size specification is referenced to 
the local body diameter ( RCORE in Item 8). The default value is 
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I 


0.025, and this value is adequate for most cases. The larger the 
core radius specified, the more the vortex effects are reduced. 
To eliminate the possibility of negating the vortex effect en- 
tirely, a maximum vortex core of .05 times the maximum body di- 
ameter is built into the code. Except in unusual situations, the 
default value is recommended, and this is the value used in the 
comparisons with data shown in a later section. 

Sample Cases 

In this section, nine sample input cases are described to 
illustrate the various program options available in VTXCLD. The 
purpose is to point out the available features and options and to 

provide a range of sample inputs which will help users prepare 

input for specific cases. Table I provides a summary of the 

sample input cases to assist the user in finding an appropriate 
example . 

Sample Case 1, Figure 10(a), is a test run based on the 3- 

caliber ogive-cylinder body in Reference 20. This run will also 
be used as the sample output case in the next section. It is 
recommended that this case be run initially to provide a check 
with the results presented in a following section. The input for 
the first sample case is described in detail, and for the fol- 
lowing cases, only the major differences are discussed. 

Input Item 1 indicates that the case consists of a body with 
a circular cross section, NCIR = 0, and a symmetric vortex wake, 
IS YM = 0. The separation criteria is laminar, NBLSEP = 1, and a 
diffuse vortex core model is to be used, NDFUS = 1. NDPHI = 1 
indicates that the two-dimensional unsteady pressure term is 
included in pressure calculation. NXFV = 7 and NFV = 0 are in- 
cluded to denote special output to be printed for seven axial 
stations which are specified in Item 30. If NFV > 0, the flow 
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field will be calculated at the 7 axial stations in addition to 
any special output. 

Input Item 2 indicates there will be four title cards input 
in Item 3. NPRNTP = 0 specifies the pressure distribution be 
output at the seven special x stations specified by NXFV , and 
NPRNTS = 1 causes a summary of the separation points to be output 
at all x stations. The flag NPRNTV = 0 results in a summary of 
the vortex field at the seven axial stations specified. NPLOTV = 
3 and NPLOTA = 1 specifies that a line printer plot of the vortex 
field is created at the seven special x stations if separation 
has occurred. NSKIN = 1 will result in an estimation of the skin 
friction on the body, and NCOMP = 0 specifies an incompressible 
calculation . 

Input Item 3 is a set of four ti tie/summary cards which are 
printed at the top of the first output page. 

Input Item 4 contains the reference information for the 
run. The reference area, REFS, is the maximum cross sectional 
area, the reference length, REFL , is the maximum diameter, and 
the moment center, XM, is located at the end of the ogive nose. 
The length of the body, SL, is 50.478 inches, and the maximum 
body diameter, SD, is 4.7 inches. 

Item 5 contains the flow conditions. The angle of attack, 
ALPHAC is 15 degrees and the Reynolds number, RE, is 440000. 
There is no roll angle for this case, PHI = 0.0. 

Item 7 specifies run parameters. The calculation procedure 
begins at x station, XI, 2.35 and continues to station, XF, 47.0 
in DX increments of 2.35 or one body radius. Separation vortices 
are initially placed a minimum 1.05 radii, EMKF, from the axis of 
the missile. The vortex reduction factor, VRF, is set to 0.6 as 
recommended for bodies which do not close themselves. 
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The only major item of interest in input Item 8 is the tol- 
erance for the vortex trajectory calculation, E5. This is set to 
0.05, a typical value. Some additional comments on this tol- 
erance are contained in the previous section on Error Messages 
and Stops. 

Input Item 9 contains the flags for the source distribution 
calculation. NSOR = -1 specifies the source distribution is 
unknown and is to be calculated by VTXCLD. NPRT = 2 provides for 
additional output regarding the source distribution. This output 
includes a comparison of the input geometry and the geometry 
calculated from the source distribution which provides a good 
means of checking the calculated source distribution. Also 
printed are the surface velocities calculated by the source dis- 
tribution. 

Input Item 12 contains NXR, the number of entries in the 
body geometry table. Items 13, 14, and 15 are the axial sta- 
tions, XR, the corresponding body radii, R, and the body slopes, 
DR, respectively. The axial stations and the radii are nondimen- 
sionalized with respect to the body length, SL. 

The final input for this case is Item 30 which identifies 
XFV, the seven axial stations where special output is printed. 

The following sample case inputs are discussed only briefly 
with some elaboration of the parameters of special interest. 

Sample Case 2, Figure 10(b), is the same as Sample Case 1 
except the source distribution is input and the angle of attack 
is changed. NSOR = 61 in Item 9 indicates that there are 61 
source locations, XSRC, and strengths, QS , input via Items 10 and 
11, respectively. The locations and strengths were obtained from 
the results of Sample Case 1. The angle of attack, ALPHAC, in 
Item 5 is changed to 20 degrees. 
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Sample Case 3, Figure 10(c), is a body with 2:1 elliptical 
cross section. This model is described in Reference 21. The 
elliptical cross section is handled by an analytical conformal 
mapping, NCIR = 1 in Item 1, and the other inputs required for 
the ellipse are the half-axis lengths, AE and BE, input in Items 
16 and 17. 

Sample Case 4, Figure 10(d), is a body of square cross sec- 
tion with rounded corners from Reference 22. A compressible Mach 
number is specified. For this cross section, a numerical con- 
formal mapping must be used, NCIR = 2 in Item 1. NCF = 0 in Item 
1 indicates that the mapping coefficients are to be calculated. 
The desired number of mapping coefficients, MNFC, is set to 20 in 
Items 19. Item 19 also specifies the number of axial stations at 
which the body is defined, MXFC = 1. MXFC = 1 is used because 
all of the axial stations are similar in shape; therefore, it is 
necessary to calculate mapping coefficients only at a single sta- 
tion. If the cross sections along the body are not similar, 
mapping coefficients must be calculated at a number of different 
sections. Item 20 is the axial station, XFC, at which the map- 
ping coefficients are calculated. In Item 21, NR = 25 is the 
number of points describing the body cross section and Items 22 
and 23 contain the horizontal and vertical coordinates of the 
cross section, respectively. 

Compressibility is included by setting the compressibility 
flag, NCOMP, to 1 in Item 1. The freestream Mach number is 
defined in Item 6, MACH =0.5. 

Sample Case 5, Figure 10(e), illustrates the required input 
for an axisymmetric body at nonzero Mach number. This case is 
Sample Case 1 with a compressibility correction. NCOMP in Item 2 
is set to 1, and the Mach number, MACH, is set to 0.8 in Item 
6. The source distribution is calculated because it changes with 
Mach number. 
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Sample Case 6, Figure 10(f), illustrates the specification 
of a transition region. The beginning and end of the transition 
region are specified in Item 7. The beginning, XTRl = 7.05, and 
the end, XTR2 = 11.75, define the region over which the separa- 
tion criteria changes from laminar to turbulent. 

Sample Case 7, Figure 10(g), illustrates the required input 
for a case in which an asymmetric wake occurs. In order to ob- 
tain an asymmetric vortex wake, the flags ISYM and NASYM in Item 
1 are set to 1. The required symmetric perturbation is specified 
in Item 8 where the initial x station for the perturbation, 
XASYMI = 0.0, and the final station, XASYMF = 5.145, are de- 
fined. The angular displacement of the separation points, DBETA , 
is 2.0 degrees over the entire region. In Item 2, NPRNTP and 
NPRNTV have been set to 1. This will cause the pressure distri- 
bution and vortex wake to be printed at each x station. 

The velocity field is calculated at 57 ( NFV) field points at 
seven axial stations, NXFV. These inputs are defined in Item 
1. The seven axial stations where the velocity field is desired 
are specified in Item 30, and the 57 field points where the ve- 
locities are calculated are included in Item 31. 

Sample Case 8, Figure 10(h), is a lobe body from Reference 
23 shown to illustrate another arbitrary cross section body. The 
noncircular input has similar format to the rounded square body 
of Sample Case 4, and the only difference is in the body cross 
section specification in Items 22 and 23 . No source distribution 
is input or calculated for this case; therefore, this situation 
approximates a two-dimensional example in which the axial coordi- 
nate represents time. 

Sample Case 9, Figure 10(i), is an example of a restart 
calculation for Sample Case 2. The case is restarted from XI - 
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21.15 and continued to XF = 47.0. Additional input required for 
this case are the overall loadings in Item 32, the vortex loca- 
tions and strengths in Items 33 and 35, and the pressure distri- 
bution in Item 37. This case is included to illustrate the means 
to restart a calculation which did not finish or one in which the 
vortex field was adjusted to include additional vortices from a 
forward lifting surface. 


Output Description 

A typical output file from VTXCLD is described in this sec- 
tion. In general, the output quantities from VTXCLD are labeled 
and each page is headed with approximate descriptive informa- 
tion. The output from Sample Case 1, Figure 10(a), requires 52 
pages, and it will be described briefly page by page. Represen- 
tative pages are shown in Figure 11. 

The first output page contains the title/summary cards, the 
reference quantities, the flow conditions, the initial run condi- 
tions, and the run flags. Page 1 is simply an echo of the input 
parameters. Also contained on this page are the force and moment 
definitions. Pages 2 and 3 echo the input body geometry, XR, R, 
and DR. The dimensional X and R values are also printed. 

Pages 4 through 10 contain information from the three-dimen- 
sional source distribution calculation. Page 4 is a summary of 
the source placement and the body geometry at these locations. 
Page 5 shows the source locations and strengths, and pages 6 and 
7 are a comparison of the input geometry shape, R/L(INPT), and 
the shape calculated from the source distribution R/L(S.D.). 
This comparison is indicative of the quality of the source dis- 
tribution. A second indication of difficulties with the source 
distribution will show up in the number of iterations required to 
find the proper shape, M-INDEX. If this number is less than four 
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or five, the chances are good the source distribution is a good 
representation of the body. If this number is large, examine the 
input and predicted shapes carefully for problems. Pages 8 and 9 
are the surface velocities calculated by the source distribution, 
and they should be examined to ensure that no exceedingly large 
velocities occur, another indication of a problem with the body 
model. Small oscillations in the velocities on the cylindrical 
portion of a body are to be expected. 

Page 10 contains a summary of source locations and 
strengths, and the values on this page are in the correct format 
for input into VTXCLD for additional cases using the same geome- 
try. See Sample Case 2 for an example. Also contained on page 
10 are the x stations at which special output is printed. 

Page 11 starts the output from the analysis part of the 
code. The first numeric line on page 11 contains the axial loca- 
tions, the local radius, the body slope, and the transition 
Reynolds number, RETR. Since axial station x = 2.35 is a special 
output station specified in Item 30, the pressure distribution is 
printed. The predicted results on pages 11 and 12 are a summary 
of the velocities and pressures on the circumference at this 
particular x station. The coordinates are dimensional, and the 
velocities are made dimensionless by the freestream velocity. 
The pressure coefficient is followed by the unsteady term, d<j>/dt. 

Following the pressure distribution is a summary of the 
section loadings and the overall loadings up to the current sec- 
tion. CN ( X ) , CY ( X ) , and CA(X) are the normal-, side-, and axial- 
force coefficients per unit length, respectively. CN, CY, and CA 
are the overall normal-, side-, and axial-force coefficients from 
the nose to the current x station, respectively. CM, CR, and CSL 
are the pitching-, yawing-, and rolling-moment coefficients, 
respectively. XCPN and XCPY are the x locations of the centers 
of pressures for the normal and side forces, respectively. 
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The next line of output is the estimation of the axial skin 
friction coefficient followed by the total axial force coeffi- 
cient made up of the pressure force and the skin friction force. 

The final block of information on page 12 contains a summary 
of the pressure distribution and the separation points at x sta- 
tion 2.35. It is noted here that no separation occurs at this 
station. 

Output pages 13 and 14 have the same format as pages 11 and 
12 so they will not be discussed in detail. The x station is 
4.7, a special output station, so the pressure distribution is 
again output. The last two lines of page 14 indicate that sepa- 
ration occurs at this axial station. The last block of informa- 
tion contains the intitial locations and strengths of the sepa- 
rated vortices. 

Page 15 shows the force and moment summary and separation 
information for x = 7.05. This format is the same as the pre- 
vious stations; however, the pressure distribution is not output 
at this station because x = 7.05 is not a special output station. 

Axial station 9.4 is again a special output station. Page 
16 is a summary of the vortex wake at this station after the 
trajectories of the individual vortices have been calculated. 
Each individual vortex in the field is described in terms of its 
location in the real plane (Y, Z) and the transformed circle 
plane (YC, ZC). The last block of information is the centroid of 
vorticity on each side of the body. This includes all the vorti- 
city, including that in the feeding sheet; therefore, the numer- 
ical centroid location may be lower than expected because of the 
effect of the vortices in the feeding sheet. 
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Page 17 has the same format as page 11 described pre- 
viously. Page 18 is a line printer plot of the cross section and 
the vortex wake. Page 19 is a summary of the separation points 
at x = 9.4, as was described for page 12. 

The remainder of the output pages, except the final three, 
repeat the formats described previously. A typical output sta- 
tion has a force and moment and separation summary like page 
15. Special output stations have a wake summary, pressure dis- 
tribution output, line printer plot, and separation summary like 
pages 16 - 19. 

The last three pages of the output shown in Figure ll(t), 
(u) and (v) give a summary of the total loadings, the vortex 
wake, and the pressure distribution. Page 50 contains a line 
printer plot of the vortex wake at the final x station and a 
summary of the total loads on the body. Page 51 is a summary of 
the vortex wake in tabular form. This information can be used 
for a restart of the calculation. Page 52 is a summary of the 
pressure distribution at the last x station; needed for a restart 
calculation . 


RESULTS 


For purposes of evaluating the accuracy and range of appli- 
cability of the subsonic vortex shedding model and associated 
computer code VTXCLD, comparisons of measured and predicted aero- 
dynamic characteristics are presented. Results for both circular 
and noncircular bodies are shown for a wide range of flow condi- 
tions including flow asymmetry and compressibility effects. 
Since the major objective is the validation of the prediction 
method, not all predicted results are in good agreement with 
experiment to illustrate to the user where problems may occur 
during general use of the code. Typical results from the predic- 
tion method follow. 
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Circular Bodies 


The prediction method was applied to an ogive-cylinder model 
(Ref. 20) in subsonic flow over a range of angles of attack to 
evaluate the pressure prediction ability of the code. The con- 
figuration has a three-caliber ogive nose and a . 7 . 7-caliber cy- 
lindrical afterbody. Circumferential pressure distributions at a 
number of axial stations are available for a range of flow condi- 
tions. 


As shown in Sample Case 1, the model is represented by 61 
point sources and sinks on the axis, and the ogive nose is used 
to close the model at the base. Because of the low Reynolds 
numbers of the experiment, laminar separation is used for all 
predictions for this model. Comparisons of measured and pre- 
dicted normal force distributions at a = 10, 15, and 20 degrees 

are presented in Figure 12. The predicted results with vortex 
wake effects included are in good agreement with experiment for 
the angles shown. The significant influence of the lee side 
vorticity is illustrated for a = 10 and 20 degrees where the 

potential flow or slender body theory results without vortex 
induced effects are shown as a dashed curve. Note that the vor- 
tex effects are small on the nose at a = 10 degrees and they 
begin to grow as the vortex increases in strength toward the end 
of the afterbody. At a = 20 degrees, the vortex effects are 

important on the nose, and they become a dominant effect on the 
cylindrical portion of the model. 

Details of the predicted characteristics are shown in Figure 
13 where comparisons of measured and predicted pressure distribu- 
tions for a = 15 and 20 degrees are presented. The predicted 

results including vortex induced effects are represented by the 
solid curves and those without vortex effects are shown as dashed 
curves. At a = 15 degrees in Figure 13(a), the vortex effects do 
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not become significant until approximately 4.5 diameters from the 
nose. Some of the roughness in the rational flow model results 
is caused by individual vortices moving too near the body surface 
during the trajectory calculation. The vortex core model does 
tend to smooth the vortex induced effects, but there can still be 
a large local effect. These local irregularities have a minimal 
effect on the integrated loads. 

The higher angle results shown in Figure 13(b) exhibit much 
larger vortex induced effects because the separation vortices are 
stronger and occur earlier on the body. The fact that there are 
more vortices with greater individual strength causes the local 
irregularities to be greater. This minor problem could be cor- 
rected to some extent by decreasing the axial interval (DX) of 
separation. The individual strengths of the vortices would de- 
crease, but there would be more vortices forming the cloud. 
Smoother pressure distributions may result at the expense of 
additional computation costs. 

The predicted vortex cloud on the ogive cylinder model at 
a = 15 degrees is shown in Figure 14. As observed above, the 
cloud does not become significant until approximately 4.5 diame- 
ters from the nose even though separation begins at one diameter 
from the nose. 

The total normal-force coefficient is shown in Figure 15 
where both potential-flow and rational-flow model results are 
compared with experiment. As noted in the pressure coefficient 
comparisons, the vortex induced effects become significant at 
a = 10 degrees and grow with increasing angle. At low angles the 
predicted results are in good agreement with the low Reynolds 
number measurements, but at high angles, the theory is in better 
agreement with the high Reynolds number results. This phenomenon 
is not completely understood at this time, and further evaluation 
with other experimental measurements is recommended. 
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The effects of compressibility on the measured and predicted 
normal force distributions are shown in Figure 16 for a = 15 
degrees. Predicted compressibility effects exhibit the correct 
trend, but they are less than those measured on the nose. The 
overshoot at the shoulder of the model is not predicted. 

Measured and predicted pressure distributions are presented 
in Figure 17. The increment in pressure coefficient due to com- 
pressibility effects appears to be in reasonable agreement with 
experiment, although there is some disagreement in the lower 
pressure regions near e > 60 degrees on the nose. This obviously 
accounts for the problem with normal force discussed above. The 
decreasing effect of compressibility along the cylindrical por- 
tion of the body is in good agreement with experiment. 

The above comparisons indicate the compressibility correc- 
tion included in the prediction method and code represents the 
measured effects; however, further verification is needed. The 
limited experimental data available are not sufficient to define 
the Mach number and angle of attack limits of the code. 

Additional comparisons of measured and predicted pressure 
distributions on other axisymmetric bodies are presented for 
further verification. For these comparisons, the 2.0- and 3.5- 
caliber ogive-cylinder models tested by Lamont (Refs. 17 and 24) 
over a wide range of angles of attack and Reynolds numbers are 
used. Only experimental data at moderate angles of attack are 
considered to guarantee flow symmetry. 

In Figure 18, measured and predicted circumferential pres- 
sure distributions on the 3.5-caliber ogive-cylinder model at 
a = 20 degrees are compared. The predictions are shown with 

(solid curve) and without (dashed curve) vortex induced effects 
included. As in the previous results, the separation-vortex ef- 
fects improve the shape of the pressure distribution curves. 
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particularly on the lee side? however, there is a general lack of 
agreement between the experiment and theory. There appears to be 
a uniform shift between the measured and predicted pressure coef- 
ficients, even on the windward side of the model where vortex 
induced effects should be negligible. 

In Figure 19, similar comparisons are shown for the 2.0- 
caliber ogive-cylinder model at a = 30 degrees. Because of the 
higher angle of attack, vortex induced effects are larger than 
before, but the character of the comparisons is nearly the same; 
that is, on the cylindrical portion of the model, the experiment 
and theory disagree by a constant increment in pressure. There 
are other areas of disagreement as it appears that the predicted 
vortex effects are often stronger than exhibited by the experi- 
ment . 


In an effort to better understand the use of the vortex 
cloud prediction method and in an attempt to resolve the differ- 
ences discussed above, the following is presented for discussion 
purposes. Since the windward meridian should have minimal vortex 
induced effects at moderate angles of attack, a comparison of 
experiment and theory along this meridian should be a measure of 
the accuracy of the linear or attached-flow part of the theory. 
In Figure 20, a number of experimental results for various ogive- 
cylinder models at a = 20 degrees are compared with a single 

predicted result. The ogive noses vary from 1.5 to 3.5 diameters 
in length (Refs. 17, 20, 24, and 25), but the pressure on the 
cylinder aft of the nose should be nearly independent of the nose 
shape. The predicted pressures are for the 3.5-caliber ogive of 
Lamont and the 3-caliber ogive of Tinling and Allen. It is ob- 
served that the theoretical results are in agreement with the 
experiments of References 20 and 25 and in disagreement with the 
experiments of References 17 and 24. The authors have no expla- 
nation for these discrepancies. 
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At higher angles of attack, the separation vorticity becomes 
asymmetric and causes both an induced normal force and side 
force. As discussed in a previous section, an asymmetric solu- 
tion develops only after the solution is perturbed, and the final 
result is somewhat dependent on the perturbation. Some prelimi- 
nary results are presented to illustrate an asymmetric result. 

Detailed flow field measurements and pressure distributions 
on a 3.0-caliber ogive-cylinder model at a = 45 degrees are 

available in Reference 26. Based on earlier experience with 
asymmetric solutions (Ref. 2), the perturbation region was lim- 
ited to the first diameter of length of the nose, and the magni- 
tude of the perturbation was varied to find an appropriate solu- 
tion. In Reference 2, the position of the shed vorticity was 
used as the criterion to evaluate the perturbation; however, in 
the present solution, a more appropriate criterion is the fre- 
guency of variation of the predicted side-force distribution. In 
Figure 21, the measured and predicted normal-force and side-force 
distributions are shown. Two different perturbations, AB = 1 and 
2 degrees, are shown, and even though the magnitude of the pre- 
dicted side force is less than that measured, the freguency of 
the side force matches best when AB = 2 degrees. This is the 

perturbation used for all the predicted results to follow. The 
predicted normal force is below the level of the measured normal 
force by a small amount, but the side force is considerably less 
than that measured. It is apparent that the vortex induced ef- 
fects are too weak for this case. 

Measured and predicted circumferential pressure distribu- 
tions at three axial stations are shown in Figure 22. Since the 
flow is asymmetric, the full circumference is shown in these 
comparisons. The major areas of disagreement are generally for 
80° < B < 160° and 200° < B < 270°; although at the aft station 
where the error in side force is greatest, the disagreement is 
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concentrated in the region 80° < 8 < 160° which causes the maxi- 
mum effect on predicted side force. The explanation of this 
problem is best accomplished by a discussion of the separation 
wake using the following figures. 

The predicted vortex cloud at x/D = 4.7 is shown, in Figure 
23 where an obvious asymmetry has developed. At this station, 
the cloud on the left side of the body has nearly broken away 
from its feeding sheet and formed a free vortex. Reverse flow 
vorticity is not included in the predicted results because of 
numerical difficulties with this particular case. The pressure 
distribution on the lee side and the reverse flow velocity field 
induced by the primary vortex field permits the calculation of 
secondary separation and the shedding of secondary vorticity. 
When these phenomena were included in the calculation, the ef- 
fects of the additional vorticity were lost when the primary 
vorticity pulled the secondary vortices into the main cloud and 
effectively eliminated their local effect on the pressure distri- 
bution. The need for the secondary vorticity is illustrated in 
the following figures. 

The measured velocity field at the same axial station is 
shown in Figure 24(a) where the individual flow vectors are 
shown. The predicted velocity vectors in the same general region 
of the lee side flow field are shown in Figure 24(b). Both mea- 
sured and predicted flow fields are similar in magnitude and flow 
directions except in the reverse flow region near the surface of 
the body. It is obvious in Figure 24(a) that a secondary vortex 
field exists in the region of interest, 120° < 8 < 170°. The 
difficulties with the predicted secondary separation vorticity 
were not resolved in the current investigation. 

The success of the previous flow field comparisons stimu- 
lated another set of comparisons to be sure the previous results 
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were not fortuitous nor accidental. Experimental flow field 
velocity components on the lee side of a 3.5-caliber ogive- 
cylinder model at a = 22.4 and 37.5 degrees are available in 

Reference 27. At the lower angle of attack, a = 22.4 degrees, 

the symmetric flow field is shown in Figure 25 for three dif- 
ferent positions in and near the vortex cloud. The location of 
the experimental results and the predicted vortex cloud are shown 
in the sketches in the figure. In each case the predicted down- 
wash velocity distribution is in very good agreement with the 
experimental results. The observed roughness when the predic- 
tions are carried out through the center of the cloud is caused 
by the influence of a discrete vortex too near the calculated 
point. 


At a = 37.5 degrees, the flow asymmetry was developed with 
a A6 = 2 degrees perturbation in accordance with the previous 

asymmetric results. The velocity field comparisons are shown in 
Figure 26, and the location of the results and the asymmetric 
vortex field is illustrated in the sketches. As above, the 
agreement is very good. The roughness of the predictions is 
again caused by the near proximity of several discrete vortices 
as seen in the sketches. 

Noncircular Bodies 

The noncircular body options in VTXCLD have been tested for 
a number of different cross section shapes, both elliptic and 
arbitrary, to verify the transformation procedures. Unfortu- 
nately, only a limited guantity of experimental data are avail- 
able for comparison purposes. A series of elliptic bodies and 
one sguare body with rounded corners are examined to illustrate 
the capabilities of the code. 
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Measured normal-force coefficients on a series of elliptic 
cross section bodies at high angles of attack are available in 
Reference 21. The cross section shapes are 1.5, 2, and 3.5:1 
ellipses with egual areas, and each model has a 3-caliber ogive 
nose and 7-caliber cylindrical afterbody. Measured and predicted 
normal-force coefficients are shown in Figure 27. Agreement 
between experiment and theory is very good for a < 20 degrees for 
all models, and the error at larger angles of attack appears to 
be nearly the same for the three cross sections. As expected, 
the nonlinear effects become more evident as the ratio of major 
to minor axes increases. 

The predicted vortex cloud patterns at the end of the body 
for each of the three cross sections are shown in Figure 28. The 
dimensionless strength of the vorticity on the positive side of 
each body is also shown in the figure. 

Though pressure measurements are not available, it is inter- 
esting to observe the change in pressure distribution along the 
elliptic bodies as the vortex strength increases. In Figure 29, 
the predicted pressure distributions around the 2:1 elliptic body 
at a = 20 degrees for three axial stations are shown. The solid 
curves from the rational flow model are compared with the poten- 
tial flow distributions, and it appears that the major nonlinear 
effects are an increase in the minimum pressure at the side of 
the body and lower pressure on the lee side of the body. 

The prediction method was next applied to a sguare body with 
rounded corners at a = 20 degrees and M b = 0.5 from Reference 
22. The model has a 3-caliber nose followed by a 10-caliber 
cylindrical afterbody. The measured and predicted normal-force 
coefficients for the body at <j> = 0 and 45 degrees roll angles are 
compared in Figure 30(a) and (b) , respectively. At 4> = 0°, the 
predicted normal force is greater than the measured value. It is 
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interesting that the measured results show very little effect of 
corner radius; that is, the smallest radius tested behaves like a 
circular cross section at subsonic Mach numbers. 

Measured and predicted results on the same body rolled 45 
degrees are shown in Figure 30(b). The agreement is very good 
for the three Mach numbers 0.5, 0.7, and 0.8. In the rolled 

condition, the effect of corner radius is significant on the 
measured normal force. This phenomena may be explained by look- 
ing at the vortex shedding results. 

In Figure 31, the predicted vortex wake on the square body 
with rounded corners is shown for <J> = 0 and 45 degrees. At 

$ = 0, separation can occur at either the upper or lower corner, 
and it may occur at both locations simultaneously if the lower 
portion of the flow reattaches on the flat side. The prediction 
method cannot handle multiple separation or reattachment; there- 
fore, it is not surprising that the normal-force coefficients in 
Figure 30(a) are in poor agreement. When the model is rolled 45 
degrees, there is only one adverse pressure region near the cor- 
ner on the side, and a single separation location is indicated. 
The prediction method appears to have little problem with this 
situation. The vortex strengths shown on each figure indicate 
that the vortex wake from the <j> = 45 degrees flow condition is 
slightly stronger than that corresponding to the unrolled condi- 
tion . 


CONCLUSIONS 

An engineering prediction method based on a rational flow 
modeling technique and the associated code VTXCLD to predict the 
vortex shedding from circular and noncircular bodies in subsonic 
flow at angles of attack and roll are described. Comparisons of 
measured and predicted aerodynamic characteristics and flow field 
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quantities are used to verify the flow model and prediction 
method for a variety of configurations under a wide range of flow 
conditions. The method has proved successful in representing the 
principal features of the complex flow field on the lee side of 
missiles at moderate angles of attack; therefore, it has applica- 
tion as an engineering or preliminary design technique directed 
at the prediction of nonlinear aerodynamic characteristics re- 
sulting from high angle of attack flows. 

The prediction method described herein has further applica- 
tion as one component of a larger prediction method for complete 
configurations consisting of an arbitrary body and multiple 
fins. The ability to model the correct flow field in the vicin- 
ity of the body leads to the capability to calculate body separa- 
tion wake induced interference effects on fins and other control 
surfaces. The vortex shedding analysis and prediction techniques 
are also applicable for use in other codes; for example, the 
methods developed in this investigation can be transferred to 
higher order codes such as panel codes. This was demonstrated in 
the supersonic vortex shedding work reported in Reference 3. 

Finally, the prediction method presented in this report must 
be considered as preliminary. Even though verification of the 
method by comparisons with experiment has demonstrated some suc- 
cess, there are configurations and flow conditions which must be 
approached with caution. For this reason, several recommenda- 
tions for improvements to the method are suggested in the follow- 
ing section. 


RECOMMENDATIONS 


In 

VTXCLD , 

u i 1. X 


the course of development and verification of the code 
several specific areas requiring additional work were 
ed. Since the additional effort was beyond the scope of 


- 81 - 


the present investigation, they are noted below as recommenda- 
tions for future work. 

The first and most important recommendation is for a thor- 
ough testing and further verification of VTXCLD. The code should 
be applied to a variety of configurations and flow conditions for 
which experimental results are available for comparison pur- 
poses. This should better define the operational limits of the 
method. The major difficulty will be the lack of useful experi- 
mental results on noncircular bodies. 

As noted in the text, use of the code at high angles of 
incidence for asymmetric vortex wake calculations presents a 
problem regarding the reguired perturbation necessary to develop 
asymmetric solutions. It is recommended that the code be applied 
to a number of asymmetric cases for which experimental data are 
available in an effort to correlate the magnitude and location of 
the perturbation in light of comparisons with data. This is the 
only method which will permit the code to be used with confidence 
for asymmetric flow conditions. 

A similar problem exists with the transition option included 
in the code. Some effort should be expended to better understand 
the effects of changing the character of separation from laminar 
to turbulent. As before, correlation of predicted results with 
experiment may provide the necessary rules for using transition 
with confidence. 

The numerical difficulty associated with secondary separa- 
tion discussed in a previous section needs further investiga- 
tion. A minor modification may correct the problem, but some 
additional effort is required to understand the source of the 
problem. 
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The prediction method may have additional uses which have 
not been explored. For example, the code has vortex shedding 
technology built into it which can be transferred to three-dimen- 
sional panel codes. The code also has potential use in providing 
separation line and vortex wake information for starting solu- 
tions for higher order prediction methods. 


Nielsen Engineering & Research, Inc. 
Mountain View, CA 94043-2287 
November 1985 
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TABLE I. SAMPLE CASES FOR PROGRAM VTXCLD 
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(b) Asymmetric separation. 


Figure 1 


Lee side vortex formation on 
an inclined body . 
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o = y + iz v=T+iA 


REAL PLANE CIRCLE PLANE 


(a) Analytical transformation procedure 





C = ior 


CL r 

= u - E J 

£=l 


£+l A, 

o <£ 


v* 


(b) Numerical transformation procedure 


Figure 2.- Conformal mapping nomenclature. 


- 90 - 



Figure 3. - Sketch of crossflow plane separation points 


i 
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Figure 4. - Body crossflow plane nomenclature. 
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Subroutine calling sequence for 
program VTXCLD. 
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Figure 6.- Subroutine cross reference list for program 
VTXCLD . 
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Figure 6.- Concluded 
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Figure 7.- Common block cross reference list for 
program VTXCLD. 
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Figure 7 . - Concluded 
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Figure 8.- Input forms for Program VTXCLD. 
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(a) Sample Case 1 

Figure 10.- Sample cases for program VTXCLD. 
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(b) Sample Case 2 
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(b) Concluded 
Figure 10.- Continued. 


- 107 - 



Item 


(1) 

1 0 

0 1 

0 1 

1 0 

8 0 

0 0 

0 0 

0 1 

(2) 

5 0 

1 0 

3 1 

0 0 

1 0 




(3) 

SAMPLE CASE 3 








NIELSEN ENGINEERING & RESEARCH. INC. 

PROGRAM VTXCLO 




3 CALIBER NOSE 



REF. 

HIT AL 138 




2 : 1 ELLIPTIC CROSS SECTION 
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INITIAL FLOW CONDITIONS ALPHAC = 15.00 DEGREES, 

PHI = 0 DEGREES 

i 

(4) 

4.4301 

2.375 
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(c) Sample Case 3 
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(c) Concluded 
Figure 10.- Continued. 
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Item 

( 1 ) 2 0 0 1 0 1 1 0 8 0 0 0 0 0 0 

(2) 5 0 1 0 3 1 0 0 1 1 

(3) SAMPLE CASE 4 

NIELSEN ENGINEERING t RESEARCH, INC. PROGRAM VTXCLD 

SQUARE WITH ROUNDED CORNERS REF. SCHNEIDER 

NUMERICAL MAPPIN6 


INITIAL FLOW CONDITIONS ALPHAC = 20.00 DEGREES, PHI = 0 DEGREES 
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(d) Sample Case 4 
Figure 10.- Continued. 
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Figure 10.- Continued- 


- 111 - 



Item 


(1) 

0 0 

0 1 

0 1 

1 0 

7 0 

0 0 

0 0 

0 

(2) 

4 0 

1 0 

3 1 

0 0 

1 1 




(3) 

SWLE CASE 5 

COMPRESSIBILITY CORRECTION 





NIELSEN ENGIhEERING & RESEARCH, INC. 

PROGRAM VTXCLD 




3 CALIBER OGIVE-CYLINDER 


REF. NASA TN D-1297 



INITIAL FLOW CONDITIONS ALPHAC = 15.00 DEGREES, 

PHI * 0 DEGREES 
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(e) Sample Case 5 
Figure 10.- Continued. 
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SAMPLE CASE 6 TRANSITION RE6I0N 
NiaSEN ENGINEERING t RESEARCH, INC. 

3 CALIBER 06IVE-CYLINDER 
INITIAL FLON CONDITIONS ALPHAC = 20.00 DEGREES, PHI = 0 DEGREES 
17.34900 4.70000 14.10000 50.47800 4.70000 

20.00000 0.00000 440000. 

2.35000 47.000 2.35000 7.05 11.75 1.05 0.0 

.04 0.0 0.0 0.0 0.0 0.0 0.0 


0 0 0 0 
0 
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-0 .92338E-04-0 .92578E-04-0 .81301E-04-0 .63716E-04-0 .46689E-04-0 .31532E-04 
-0 .20487E-04-0 .12619E-04-0 .76726E-05-0 .43306E-05-0 .27866E-05-0 .10763E-05 
-0.14198E-05 
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(f) Sample Case 6 


Figure 10.- Continued. 
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(f) Concluded 


Figure 10.- Continued. 
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*mE CASE 7 ASYMMETRIC SEPARATION 
NIELSEN ENGINEERING t RESEARCH, INC. PROGWM VTXCLD 

3 CALIBER OGIVE-CYLINDER REF. NSHC TR82-384 

ASYMMETRIC SEPMATION CASE 
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(g) Sample Case 7 
Figure 10.- Continued. 
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(g) Concluded 
Figure 10.- Continued. 
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SAMPLE CASE 8 
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LOSE BODY 
NUMERICAL HAPPING 


INITIAL 

FL0H CONDITIONS ALPHAC = 15.00 
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0.08840E+00 0.10026E+00 0.11049E+00 0.11882E+00 0.12508E+00 0.12926E+00 
0.13147E+00 0.13194E+00 0.13097E+00 0.12891E+00 0.12607E+00 0.11400E+00 
0.10000E+00 0.08365E+00 0.07847E+00 0.07212E+00 0.06431E+00 0.05476E+00 
0.04333E+00 0.03011E+00 0.01546E+00 O.OOOOQE+OO 
-0 .17874E+00-0.17809E+00-0 .17612E+00-0 .17275E+00-0 .16786E+00-0 .16134E+00 
-0 .15312E+00-0 .1431 9E+00-0. 131 68E+00-0 .11882E+00-0 .10496E+00-0 .09051E+00 
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(h) Sample Case 8 


Figure 10.- Continued. 
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SIMPLE CASE 9 

NIELSEN EN6IICERIN6 & RESEARCH, INC. 
3 CALIBER OGIVE-CYLINDER 
RESTART CALCULATION 


0 4 0 9 0 9 

0 1 0 

PROGRAM VTXCLD 
REF. MSA TN D— 1297 
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FLON CONDITIONS ALPHAC = 20.00 

DEGREES, PHI = 0 DEGREES 
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0.48256E+00 0.51747E+00 0.55238E+00 0.58730E+00 0.62221E+00 0.65712E+00 
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(i) Sample Case 9 
Figure 10.- Concluded. 
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Figure 11.- Continued 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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Figure 11.- Continued 
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Figure 11.- Continued 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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FORCE AND MOMENT COEFFICIENTS - PRESSURE INTEGRATION 
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Figure 11.- Continued. 
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Figure 11.- Continued 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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Figure 11.- Concluded. 













G 

O 

•H 

■P 

G 

;0 

-P 

cn 

*H 


Q) 

G 

G 

in 

in 

<D 

U 


LD 
i— I 

II 

a 

rtf 


0) 

-p 

o 

■H 

TJ 

CD 

G 

04 

G 

fd 

TJ 

CD 

G 

G 

tn 

<d 

CD 

s 


G 

(d 

G 

O 


*■ 

CO 


CD 

G 

G 

•r+ 

p4 


144 


ogive-cylinder model. 
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Figure 15.- Measured and predicted normal-force 
coefficient on an ogive-cylinder 
model. 
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Experiment, Ref. 20 
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Figure 20.,- Measured and predicted pressure 

distribution on the windward meridian 
of several ogive-cylinder models at 
a = 20°. 
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Comparison of measured and predicte< 
distributions of normal force and s 
force on an ogive-cylinder model at 
= 45 ° . 
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Figure 23,— Predicted asymmetric wake 
at x/D — 4.7 on an ogive-? 
cylinder model at a = 45° 
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Figure 24.- Measured and predicted velocity field at x/D 
on an ogive-cylinder model at a = 45°. 
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Figure 28.- Predicted vortex cloud wakes on a 

series of elliptic cross section bodies 
at x/D = 10, a = 20°. 


- 160 - 




- 161 - 


















1. Report No. 

NASA CR-4037 


2. Government Accession No. 


3. Recipient's Catalog No. 


5. Report Date 

January 1987 


6. Performing Organization Code 

2021 /C 


8. Performing Organization Report No. 

NEAR TR 354 


10. Work Unit No. 


11. Contract or Grant No. 

NASl-17962 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 

533-02-01-04 


4. Title and Subtitle 

Prediction of Vortex Shedding From Circular and 
Noncircular Bodies in Subsonic Flow 


7. Author(s) 

Michael R. Mendenhall and 

Daniel J. Lesieutre ' 


9. Performing Organization Name and Address 

Nielsen Engineering and Research, Inc. 

510 Clyde Avenue 
Mountain View, CA 94043 

1 2. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


15. Supplementary Notes 

NASA Langley Research Center 
Technical Monitor: Robert M. Hall 

16. Abstract 

An engineering prediction method and associated computer code VTXCLD are 
presented which predict nose vortex shedding from circular and noncircular 
bodies in subsonic flow at angles of attack and roll. The axisymmetric body 
is represented by point sources and doublets, and noncircular cross sections 
are transformed to a circle by either analytical or numerical conformal trans- 
formations. The leeward vortices are modeled by discrete vortices in cross- 
flow planes along the body; thus, the three-dimensional steady flow problem is 
reduced to a two-dimensional, unsteady, separated flow problem for solution. 
Comparison of measured and predicted surface pressure distributions, flow- 
field surveys, and aerodynamic characteristics are presented for bodies with 
circular and noncircular cross sectional shapes. 


17. Key Words (Suggested by Author(s) ) 

Vortex shedding 
Circular bodies 
Noncircular bodies 
Discrete vortices 

19. Security Oassif. (of this report) 

Unclassified 


18. Distribution Statement 


Unclassified - Unlimited 


Subject Category 34 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages 
168 


22. Price 


N-305 


For sale by the National Technical Information Service, Springfield. Virginia 22161 











